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Abstract £g 

Solutions to the nuclear many-body problem rely on effective interactions, and in general ef- 
fective operators, to take into account effects not included in calculations. These include effects ^ 
due to the truncation to finite model spaces where a numerical calculation is tractable, as well as p£ 
physical terms not included in the description in the first place. In the no-core shell model (NCSM) *-p 
framework, we discuss two approaches to the effective interactions based on (i) unitary transfer- ^ 
mations and (ii) effective field theory (EFT) principles. Starting from a given Hamiltonian, the 
unitary transformation approach is designed to take into account effects induced by the truncation 
to finite model spaces in which a numerical calculation is performed. This approach was widely 
applied to the description of nuclear properties of light nuclei; we review the theory and present 
representative results. In the EFT approach, a Hamiltonian is always constructed in a truncated 
model space according to the symmetries of the underlying theory, making use of power counting to 
limit the number of interactions included in the calculations. Hence, physical terms not explicitly 
included in the calculation are treated on the same footing with the truncation to a finite model 
space. In this approach, we review results for both nuclear and trapped atomic systems, for which 
the effective theories are formally similar, albeit describing different underlying physics. Finally, 
the application of the EFT method of constructing effective interactions to Gamow shell model is 
briefly discussed. 



1 Introduction 

It is generally accepted that the nuclear problem is pathologically complicated. On one hand, the 
inter-nucleon interaction is non-central, non-local and unconstrained at short distances, and, on the 
other hand, the complexity of solving the many-body Schrodinger equation increases dramatically with 
the number of nucleons. Phenomenological [1, 2] and one-boson exchange models [3] of the inter- 
nucleon interactions have provided guidance and have proved successful in applications to light nuclei. 
However, a deeper understanding of the interactions between nucleons has been achieved using effective 
field theories (EFTs) [4, 5], which provide interactions consistent with the symmetries of the underlying 
theory of the strong interactions, QCD. EFT has several advantages: eliminates the model dependence 
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(e.g., there are an infinite number of nucleon-nucleon interactions which are phase-shift equivalent, 
but many-body solutions depend upon the particular choice), explains naturally the hierarchy of the 
nuclear interactions, provides a framework for error estimation, and can mitigate shortcomings with the 
description of low-momentum observables, such as multipole transitions. 

For light nuclei, several methods are now available for solving the nuclear few-body problem, such as 
Green's function Monte Carlo (GFMC) [6, 7], Fadeev [8], Fadeev-Yakubovsky [9], (effective-interaction) 
Hyperspherical Harmonics [10, 11], and the no-core shell model (NCSM) [12, 13], and they all agree 
within the statistical or numerical errors inherent to each method [14]. Coupled-cluster (CC) expansion 
has been applied to the description of (nearly) closed-shell nuclei, from 4 He to 40 Ca [15, 16]. In most 
cases, the theoretical description of the experimental data is excellent. Among the few-body methods, 
while subject to greater errors, the NCSM, CC, and the auxiliarly-field difusion quantum Monte Carlo 
(AFDMC) [17], are the only methods with the potential to be extended to heavier nuclei with fewer 
restrictions. Thus, NCSM can handle both local and non-local interactions unlike GFMC that has 
difficulties in the presence of strong non-localities. AFDMC has been restricted to more schematic 
interactions like Argonne v$, and to date only an approximate inclusion of three-body forces has been 
possible in CC. It is also conceivable that a lattice approach like the one in Ref. [18] could be applied 
to heavy N = Z nuclei. 

NCSM is a direct numerical diagonalization approach, in which the Schrodinger equation is solved by 
expanding the system wave function in a many-body basis. The many-body basis is truncated to a finite 
size, so that a numerical diagonalization can be performed. This is similar to the phenomenological 
shell model, shortly described in Sec. 2, which deals with even more drastic truncations of the available 
model spaces. However, in contrast with the phenomenological version, in NCSM one starts with a given 
Hamiltonian (which can contain two- and three-body interactions) defined in an infinite Hilbert space 
and one derives an effective interaction in a given model space via a unitary transformation. We review 
the traditional derivation for effective interactions, based on the Okubo-Lee-Suzuki transformation 
[19, 20, 21, 22], as well as effective operators, and present selected representative results in Sec. 3. 

Applications of few-body methods like GFMC and NCSM has allowed tremendous progress in the 
last decade in our understanding of how nuclear structure arises from the properties of the interactions 
among nucleons inside a nucleus. Often, these interactions are modeled only in terms of ad hoc po- 
tentials, which are not necessarily rooted in QCD. EFT provides the framework to construct nuclear 
potentials that respect the symmetries of QCD and produce observables in a systematic and controlled 
expansion in a small parameter [4, 5]. In such an approach, the shape of the nuclear interactions, 
including the range of pion exchange, is restricted and the complicated short-range physics is encoded 
into a number low-energy constants, which, so far have been fitted to data. However, given the advances 
in lattice QCD simulations [23], it will be possible in the near future to directly compute these from 
the underlying theory. 

While the EFTs can shed light on the properties of the inter-nucleon interactions, the solutions 
to many-body problems, especially those involving a large number of nucleons, rely on the increasing 
computational power, supported by the development of sophisticated numerical algorithms. In general, 
the development of many-body methods are considered independent from the derivation of the interac- 
tions. To some extent, this is a justified approach especially if one needs to asses the reliability of the 
many-body method, so that in conventional approaches the two issues have been considered decoupled: 
the effective interaction has been seen as merely an input to the many-body codes, with minimal mixing 
between the two worlds, required for the determination of the three-nucleon interaction parameters. 

With the development of EFT, it has been recognized that one has never access to the "full" Hilbert 
space associated with the quantum system. It becomes therefore necessary to truncate the Hilbert 
space in order to exclude those states that could not be constrained experimentally or theoretically. In 
other words, the interactions are only defined in the context of a model space, and often independent 
of whether the excluded space physics is known or not. Much in the same way one is forced to restrict 
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the calculation to a tractable model space in NCSM. This has motivated the development of a hybrid 
approach, in which the interactions are derived using the principles of EFT directly in NCSM model 
spaces. In this paper, we review the efforts to combine the power of EFTs with the NCSM, with the 
goal of not only extending the benefits of QCD compatible solutions to a larger class of nuclei, but also 
of eliminating the model dependence and mitigating some shortcomings of the conventional NCSM. 
Reference [24] was a first attempt at generating effective interactions in NCSM restricted spaces, by 
fitting two- and three-body parameters directly to levels in light nuclei. An alternate method of fitting 
the leading order (LO) low-energy constants for the two-body interaction was presented in Ref. [25], 
where a physical trap was added. The method has been later extended in order to include corrections 
beyond leading order, with applications to few-body systems of cold atoms [26, 27] and nucleons [28]. 
We present in detail the two methods in Sec. 4. Furthermore, the approach is general and has been 
recently applied to the derivation of effective interactions used in Gamow shell model calculations, as 
illustrated in Sec. 5. We conclude and discuss further applications in Sec. 6. 

2 Phenomenological shell model 

Although the main interest of the present paper is NCSM, we find useful to briefly review the con- 
ventional shell model, which is the phenomenological method of choice for the description of medium 
nuclei. To some extent, the theoretical derivation of the effective interaction is even more challenging 
than in the NCSM, because a much larger amount of correlations have to be included in the effective 
Hamiltonian. 

In the phenomenological shell model, all the correlations between a reduced number of nucleons are 
included. The tradeoff is that most of the nucleons are considered inert. Thus, the main assumption is 
that the wave function factorizes in a component describing a closed shell nucleus, and a few valence 
nucleons: 



and valence nucleon wavefunctions respectively. A necessary condition for the factorization (1) to work 
is that the core subcomponent of the system have no low-lying excited states. While this is condition 
is satisfied for p-shell nuclei, because the core nucleus 4 He has no low-lying excited states below 20 
MeV, the condition could be badly violated for other cases. For example, 16 0, the core for so-called 
sd-shell nuclei (mass number between 20 and 40), has low-lying excited states around 4 MeV. In such 
cases, in order to keep the simple factorization of the wave function (1), one needs to include in the 
effective operators the effects of the nucleons being excited out of the closed-shell core, process called 
core polarization [29]. Assuming that one can separate the nuclear Hamiltonian in core and valence 
contributions 



the eigenvalue equations decouple in a part that involves only the core (and is assumed solved) 





(2) 




(3) 



H val \* wd ) = E val \* vai ), 



(4) 



which constitutes the main object of investigation in the phenomenological shell model. 
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Even if we assume the wave-function factorization produces small enough errors, one needs further 
simplifications since the valence particle wavefunction \ty val ) given by Eq. (4) allows particles to scat- 
ter into an infinite number of states above the core. In order to make calculations tractable, in the 
phenomenological approach one considers that the valence space spans only a small enough number of 
states. For example, for p-shell nuclei, a good description of a large number of low-lying states can be 
obtained if the lowest p-shell single particle states are accessible to the valence nucleons. The tradeoff 
is that only a subset of low-lying states can be described in such a restricted space, even if all miss- 
ing correlations are included in the effective operators. Thus, for even- A nuclei one can only describe 
positive-parity states, while for odd- A nuclei only negative-parity states are accessible. In order to 
describe the remaining states, one has to work in valence spaces that include multiple shells. 

In order to derive the interaction in the restricted model space, we consider two interacting particles. 
Let |r) be the wave function so that the Schrodinger equation is 

(H q + V)\t) = E t \t), (5) 

where H is the kinetic energy term, and V can be, in principle, the bare nucleon-nucleon interaction. 
The last equation transforms into two equations if we multiply on the left by P and Q, the projectors 
inside and outside the valence space respectively (by definition, the projector P projects into the valence 
space, i.e., P\t) = \^ val )): 

{-E T + H + PVP)P\r) = -PVQ\t), (6) 

(-E T + H + QVQ)Q\t) = -QVP\r). (7) 

Solving the latter for the wave function outside the model space Q\t), and plugging into the former, 
one obtains a Schrodinger equation inside the valence space [30]: 

(H + PV e f f P)P\r) = E T P\r), (8) 

with the effective nucleon-nucleon interaction in the valence space given by 

V$ = V + VQ Er _ H l_ QVQ QV. (9) 

This can be rewritten as an integral equation for [30]; furthermore, in the particular case when the 
Q operator corresponds to the two-particle excited states outside the model space (which brings the 
largest contribution), one can show that the equation for V e ff is similar to the equation for Briieckner 
G-matrix [31]. Therefore, computation of effective interaction reduces to calculation of the G-matrix. 

In general, the effective interaction obtained through the procedure outlined above does not provide 
good description of the low-lying states. Therefore, one further adjusts the two-body interaction using a 
fitting procedure so that one obtains a correct description of ground and excited states for a large number 
of nuclei. Very often, even if one starts with a realistic two-body potential, for a better description of 
the experimental data, one can completely ignore the analytical form of the nucleon-nucleon potential 
V(fi,fj). Thus, the effective TBME approximation consists in assuming that each matrix element 
is assumed a parameter that can be adjusted in order to obtain the experimental spectra of nuclei, 
and the G- matrix is usually the first step of the iteration. Because it is analytically convenient, one 
usually chooses single-particle states corresponding to the harmonic oscillator. This has no relevance 
for diagonalization, but it becomes important for calculating transition strengths when one has to 
evaluate matrix elements of different operators, as the reduced matrix elements of the corresponding 
operator depend upon the radial structure of the single-particle wave function. In general, however, in 
the phenomenological shell model, no diagrammatic basis for generating the effective operator exists. 
Empirical renormalization of transition operators must also be introduced, obscuring the underlying 
physics and undercutting the phenomenological shell model as a predictive tool. 
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However, the phenomenological shell model is considered a success. In the full sd space the parame- 
ters are the three single-particle energies (the spin-orbit interaction removes the degeneracy) £i Sl/2 , £od i/2 
and e 0ds/2 plus 63 TBME. The Wildenthal "USD" interaction [32] is fitted to reproduce 447 ground and 
excited states of sd shell nuclei (A = 17 — 40). In the space spanned by the lpi/2 — 1^3/2 — O/5/2 — O/7/2 
single particle states outside a 40 Ca inert core, the monopole-modified Kuo-Brown "KB3" interaction 
[33] is very successful in reproducing experimental features; this interaction was derived starting from 
the G-matrix as first approximation, and then fitting the matrix elements to describe experimental data. 
Further refinement of the existing interactions to include more experimental information in the fitting 
procedure have allowed even more precise description of nuclear properties both in the sd [34] and pf 
[35] shells. A more systematic approach to the derivation of effective interactions [36] and operators 
[37] in one single shell has been implemented within the NCSM framework, and will be discussed in 



Phenomenological shell-model inspired methods have been applied to the description of trapped 
gases, even though the model spaces are not so drastically restricted [38]. Given the fact that the off-shell 
properties of the NN interactions cannot be constrained experimentally, Johnson has proposed a unitary 
transformation that preserves the on-shell properties, but is fitted, similar to the phenomenological 
interactions, to give the correct binding energies for select few-body systems, thus minimizing the 
expectation value of the induced few-body forces. Using this method, several few-body systems not 
included in the fit have been predicted with good accuracy in small model spaces [38]. Thus, it is 
possible that such an approach could be used for the derivation of effective interactions in small model 
spaces for heavy nuclei, although its predictability should be properly investigated. 



NCSM is a few- and many-body technique specialized in solving the Schrodinger equation by direct 
diagonalization in a restricted basis set constructed with harmonic oscillator (HO) wave functions. 
The method has its roots in the phenomenological shell model described in the previous section. The 
advantage of this approach, however, is that the results are independent upon any parameters intrinsic to 
the calculation, such as the HO frequency of the single particle states used to construct the many-body 
basis, even though in practice this is achieved only for small systems. 

In NCSM all the nucleons are allowed to interact. Formally, the non-relativistic intrinsic Hamiltonian 
describing the system of A protons and neutrons writes 



where is the NN interaction, which depends only on the relative coordinates between the i and j 
particle, V^ N the NNN interaction, pi the momentum of particle i, and Mm the nucleon mass. Usually, 
V NN is fitted with high accuracy low-energy observables (phase-shifts and deuteron properties), while 
the NNN interactions are adjusted to reproduce properties of the three-body system (tritium, nucleon- 
deuteron scattering). This Hamiltonian can be numerically diagonalized in an appropriate basis. 

There are two equivalent approaches to solving the many-body problem, depending on how the basis 
is set up. In the first approach, one considers the relative (or Jacobi) coordinates, defined, for example, 
in terms of differences between the CM positions of sub-clusters within the A-hody system: 



Sec. 3.5. 



3 No-core shell model 




(10) 
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(ri + r 2 H h r*A-i) - 



(11) 



Subsequently, the internal coordinates are described with HO wave functions [39]. For example, the 
basis states for the three-nucleon system are 



which have the spatial part constructed using HO wavefunctions in ^ and £ 2 with quantum numbers 
ni, h and n 2 , / 2 , respectively, with the angular momentum coupled to L, while the spin (isospin) part 
is constructed by coupling first two spins (isospins) s — 1/2 (t — 1/2) into spin (isospin) s 2 (t 2 ) and 
then a third spin (isospin) s = 1/2 (t = 1/2) to total spin S (isospin T). In Eq. (12), A stands for 
the operator that antisymmetrizes the three-particle wavefunction. Details on the construction of a 
fully antisymmetrized basis can be found in Ref. [39], which include generalization to more than three 
particles. Because in this method only relative coordinates are involved, the solutions are translationally 
invariant. However, due to the complexity of the anti-symmetrization, this approach is not efficient for 
nuclei with A > 4. 

For heavier nuclei, it is more efficient to work with a Slater determinant (SD) basis, in which 
antisymmetric many-body states are constructed from single HO basis states. In such an approach, 
center-of-mass (CM) excitations can mix with the internal degrees of freedom. However, if one performs 
a truncation on the number of energy quanta above the non-interacting minimum configuration, the 
CM wave function factorizes exactly, so that the results are free of spurious modes. In particular, by 
adding 



to the Hamiltonian (10), with u> the HO frequency and j3 a numerical parameter of order 10, one can 
ensure that for the low- lying states, the CM is in the lowest HO. Hence, for operators depending only 
upon the intrinsic coordinates, the matrix elements are free of CM contributions. While in this case 
the anti-symmetrization is trivial, the size of the many-body model space increases very quickly with 
the number of particles and the number of single-particle states. Even so, the diagonalization in a SD 
basis remains the method of choice for A > 4, and a lot of effort is dedicated toward the development 
of algorithms to handle larger and larger dimensions. 

In principle, NCSM provides exact solutions to the Hamiltonian in Eq. (10). In either the trans- 
lationally-invariant or the SD approaches, for a complete calculation, an infinite number of HO states 
should be included in calculations, and therefore the dimension of the many-body basis is infinite. Since 
the NCSM involves a numerical diagonalization, the many-body basis has to be truncated, and hence a 
method to include correlations left out by truncation is required. In NCSM, the truncation is determined 
by N max , which is the number of u excitations above the minimum non-interacting configuration. Thus, 
the model space is defined by N max and the frequency u>. Note that if the truncation in the SD basis, 
given by N max , is related to the truncation in relative coordinates, so that both include the same number 
of excitations on top of the minimum configuration, calculations by the two means produce the same 
results for the same truncations. 

The intrinsic properties of the A-body system are not affected by the addition of a Hamiltonian 
term that depends only on the CM coordinates. Thus, following Lipkin's idea [40], we add a HO CM 
Hamiltonian, so that the new Hamiltonian reads 




(12) 
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H A = H mt + HS% = E^+E V - ^ + E + (13) 

i=l i>j'=l ^ ' i>i>fe=l 

where hi stands for a single-particle HO Hamiltonian of frequency u. As one subtracts the CM term 
in the final many-body calculation, it does not introduce any net influence on the converged intrinsic 
properties of the many-body calculation. Furthermore, this addition and subtraction does not affect 
our exact treatment of the CM motion. However, this procedure is an essential step in the derivation 
of the effective interaction within the NCSM approach. 



3.1 Unitary transformation method for effective interactions 

In most of the approaches to solving the nuclear many-body problem, the inter-nucleon interactions are 
considered input. Early applications of NCSM have used as input phenomenological NN interactions, 
such as the local Argonne vis [1, 2] or the non-local charge-dependent (CD) Bonn interaction [3]. In 
such interactions, the short-range part has been parametrized as a repulsive code of the order of 1 
GeV. As a consequence, even low-lying states have considerable high-momentum components, making 
the convergence in finite model spaces very difficult. The unitary transformation approach takes into 
account the eliminated states to produce an effective interaction for the truncated model space. Any 
truncation induces up to A-bodj interactions, so an effective interaction that reproduces exactly the 
low-lying spectrum is as difficult to derive as solving the A-bodj problem. Hence, the so called cluster 
approximation has been derived, where the interaction is obtained for a < A particles, and the result 
used in the A-body problem. However, in finite models spaces, the addition of the CM procedure 
discussed at the end of the previous section introduces a pseudo-dependence upon the HO frequency 
ou, and the cluster approximation described below will sense this dependence. In the largest model 
spaces, important observables manifest a relatively weak dependence of the frequency u and the model 
space size. In the following, as it is customary in NCSM applications, we make the distinction usual 
distinction in the traditional NCSM between the "bare" interaction, which is the original interaction 
whose full space full Hilbert space solution we seek, and the "effective" interaction, which is derived 
in a small model space. We will see that in the EFT inspired approach discussed later that such a 
distinction is artificial. 

Perhaps the most general derivation of the unitary transformation approach was presented in Refs. 
[41, 42], which we closely follow in this paper. 

As discussed previously, for a given Hamiltonian H int , in particular given by Eq. (10), the goal is 
to solve the Schrodinger equation 

H int \$) = E*\Q) (14) 

in a finite model space, by direct numerical diagonalization. Hence, in the NCSM approach one divides 
the full Hilbert space into a model space, with associated projection operator P (P 2 = P), and a 
complementary, excluded space, with the associated projection operator Q (Q 2 — Q, P + Q — 1, and 
PQ = QP = 0). In principle, the excluded space Q is infinite, but in practical realizations is taken 
to be as large as possible and still be numerically tractable. In addition, we introduce a similarity 
transformation operator X (not necessary unitary) with the goal to perform many-body calculations in 
the model space P, using a transformed Hamiltonian H, 

n = XH mt X~\ (15) 

so that a finite subset of eigenvalues of the initial Hamiltonian Hi nt in Eq. (10) are reproduced. This is 
a general approach, which can be applied to non-Hermitian Hamiltonian operators that can arise, for 
example, in the context of boson mappings. 
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In order to derive the conditions to be imposed on the transformation operator X, we revisit the 
Feshbach projection formalism applied to the Schrodinger equation 

H\V) = Ey\V). (16) 

We note that, for the same energy eigenvalue, the wavefunctions |$) and are connected via the 
unitary transformation 

|tf)=X|$). (17) 

In general for non-Hermitian Hamiltonians the left and right eigenvectors are not related simply by a 
Hermitian conjugation. Nevertheless, we have the freedom to choose a normalization so that {^e\^e) — 
1, where (&e\ is the left eigenvector corresponding to the eigenvalue Ey. From Eq. (16), we obtain 
immediately two coupled equations for the components of the wave function in the models space P, 

PHP\^) + PHQ\^) = E^P\^), (18) 

and in the excluded space 

QrHP\*) + QHQ\*) = E*Q\V). (19) 
Solving formally for Q\^>) in Eq. (19) 

QW = E^QHQ QHP ^ (20 > 

and inserting the result into Eq. (18), one obtains immediately that the effective Hamiltonian in the 
model space can be expressed as 

H eff = pup + puq e ^ _ X QnQ Qnp, (21) 

which is manifestly energy and state dependent. While approaches in which the energy dependence 
have been proposed, but their application has been limited to two- and three-particle systems [43, 44]. 
Using Eq. (21), we can eliminate at once the state dependence and energy dependence if we impose 
one of the following decoupling conditions 

QUP = 0, (22) 

or 

PUQ = 0. (23) 

We note, however, that the former condition also ensures that the Q-space component of the wave 
function vanishes, although this is not true for its complementary left eigenstate. Moreover, as it 
will become clear in the derivation of the effective operators below, both conditions have to be satisfied 
so that one obtains energy-independent effective operators corresponding to other observables besides 
the Hamiltonian. 

In a consistent approach general operators O are transformed by the same transformation operator 
X, O = XOX~ l . In this case, one needs to compute a matrix element of the form (<E>|(9|\I f ), where (<E>| 
corresponds possibly to another left eigenvector of H. Similar to Eq. (20), the Q-component of the left 
eigenstate (<3>| can be written 

^ Q= w pHQ E^om- (24) 

one can immediately extract the expression for the effective operator in the model space P 

°>» = POP + PnQ E^QHQ Q ° P + P ° Q E^QHQ QnP 

+PUQ E^mi QOQ E^QHQ QnP < 25 > 



8 



As advertised, in order to obtain an energy-independent expression for a general effective operator one 
needs to construct the transformation operator X so that both decoupling conditions (22) and (23) 
are satisfied. Consequently, both left and right P eigenstates of the transformed Hamiltonian % have 
components only in the model space. Additional subtleties exist within this effective operator approach 
[45]. 

In order to obtain the transformation operator X, let us consider first an operator Cj that satisfies 
the following condition: 

Cj = QCjP. (26) 
In terms of Cj, let us further assume that the transformation operator X can be written as 

X = exp(-w) = P + Q-Cj, (27) 

so that the decoupling condition (23) is satisfied automatically. The remaining decoupling condition 
(22) provides a quadratic equation that determines the auxiliary operator Cj: 

QHP = QHP - QCjHP + QHCjP - CjHCj = 0. (28) 

Finally, we obtain at once the expressions for the effective Hamiltonian in the model space 

Heff = PHP + PHCjP, (29) 

and similarly for the effective operator 

O eff = POP + POCjP. (30) 

Manifestly, the Hamiltonian (29) is non-Hermitian, even if the original Hamiltonian is Hermitian, and 
the effective operators given by Eq. (30) also change symmetry properties under a Hermitian conjuga- 
tion. 

As noted before, it is desirable to obtain Hermitian Hamiltonians in the model space, as they are 
easier to handle numerically. Thus, assuming that X = exp(G), where G is anti-Hermitian (G^ = —G), 
one obtains a considerably more complicated equation that determines G (see Eq. (2.19) in Ref. [22]). 
However, one can show that one can derive a relation between G and Cj [22] 

oo 1 

G = arctanh(w - Cj ] ) = ^(-l) n w(tfw) n - h.c, (31) 

where we have used the property CC k = 0, for k > 1. Using G given by the previous equation, we obtain 
the unitary transformation that ensures the decoupling of the model and excluded spaces 

X = (l+u-tf)(l + utf + tfu). (32) 

If we calculate QHP we obtain for Cj^ an equation which is exactly the conjugate of (28). Finally, 
because Cj^Cj and CjCj^ act only in the model space P and in the excluded space Q respectively, the 
expression for the effective Hamiltonian in the model space is 

n eff = P + " f H P + " , (33) 
St VP + Cj^Cj VP + Cj^Cj 

while the expression for general effective operators is similar 

O eff = P + ^ O P + " ■ (34) 
;/ VP + CMj VP + CM 
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The construction method we have presented here is not unique. In Ref. [41], the derivation is a bit 
more involved, using an additional operator. However, the end result is the same, with the effective 
operators in the model space given by Eqs. (33) and (34). We should, however, point out that our 
derivation does not exclude the existence of a different unitary transformation that would ensure the 
same decoupling between the model and excluded space. 

The fact that one usually chooses to work with the Hermitian Hamiltonian (33) is a matter of con- 
venience. As discusses above, there exists at least one similarity transformation, Eq. (27), that ensures 
the decoupling condition, and there could be more. Nevertheless, as long as the transformations are 
computed exactly, they all produce the same spectrum and properties of the system. The wave functions 
will be different, but this is not alarming since they are not observables. As we will see in the following, 
however, one always uses an approximation, so that the decoupling condition is only approximate, and 
the exact decoupling is achieved only in a large model space limit. It is therefore conceivable that 
one transformation could be more suitable than another in the sense of faster convergence to "full" 
Hilbert space results. Where relevant, we will discuss further implications of the particular choice of 
transformation. Hitherto there is no approximation in our derivation. 

A great deal of effort has been directed toward the calculation of the operator Co [21, 22, 46, 47, 48], as 
it is an essential step in the derivation of the model space effective operators. Two iterative solutions have 
been devised: one that converges to the states with the largest P-space components and is equivalent 
to the solution of Krenciglowa and Kuo [46], and another which converges to states lying closest to a 
chosen parameter appearing in the iteration procedure [21, 22]. A more efficient method was introduced 
later by Navratil and Barrett [47, 48] . It relies on the fact that the components of the exact eigenvectors 
in the complementary space are mapped into the model space. Thus, going back Eq. (17) and using 
the definition of the similarity transformation operator (27) we obtain 

|tt) = P|$)+w|$), (35) 

taking also into account the decoupling condition (22). Equation (35) represents the formal proof that 
Co maps back from the model space the excluded component of the wave function. Hence, if we choose 
a set /C of exact eigenvectors of the Hamiltonian in the full space, one obtains immediately that the 
matrix elements of the operator Co are given by [12, 47, 48] 

(a Q \Co\a P ) = J2( a Q\^k)(^k\a P y\ (36) 

k&K 

where \ap) and \o.q) are the basis states of the P and Q spaces, respectively. The dimension of the 
subspace JC is equal with the dimension of the model space P and we assume that the overlap matrix 
between the chosen eigenvectors and the basis states in the model space is non singular. In the next 
subsection, we will present a practical implementation of Eq. (36). 

3.2 The cluster approximation 

A closer examination of Eq. (36) reveals that the computation of the matrix elements of the operator 
Co requires the knowledge of a set of A-bodj eigenvectors, which is in itself the goal of calculation. 
Moreover, it shows that the effective Hamiltonian is a A-body operator. Hence, an exact calculation of 
the unitary transformation is as difficult as obtaining the many-body solution. The practical solution 
is the so-called cluster approximation, in which the unitary transformation operator is calculated for 
a < A particles, and an effective a-body interaction in the model space is used to solve the full A-body 
problem. Only the two- and three-body cluster approximation has been implemented up to now. 

Let us consider for the moment that one can neglect for now three- and higher-body forces and let 
us consider the a-body cluster approximation. In this case, from the A-body Hamiltonian in Eq. (13) 
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one considers the a-body problem 



h^ = ±h l+ ± (V, 

i=l i>j=l ^ 



M N cu 



2 



n 



f,) 2 . (37) 



The a-body Hamiltonian h a is further separated into a relative contribution and a pair CM contribution, 
so that the previous equation can be cast as 

= h% + tff M , (38) 
with the CM of the a-body system decoupling from the relative Hamiltonian given by 



a 

i<j=l 



2Mm 2 V A 



+ T j V ij . (39) 



i<j 



We assume that Eq. (39) can be solved in a space large enough space that could be considered the full 
Hilbert space (although in practice is always truncated), obtaining a set of eigenvectors |^^), which 
in turn determines by means of Eq. (36) the operator co^ a \ Hence, the effective a-body interaction in 
the model space will be given by the analogous of Eq. (33) 

y(q) = P ° + " (a)t P + " (a) f ht (40) 

y/P a + U^W^ ^P a + U}( a W a ) ~t 

so that the effective Hamiltonian for the A-body problem becomes 

{ a ) { 2 ) " < "-'- 1 

There is no summation over a. Note that the decoupling conditions (22) and (23) are now valid in the 
a-body space and not the full A-body space. 

In particular, in the lowest approximation, i.e., the two-body cluster approximation, the relative 
two-body Hamiltonian (39) writes 

ft "< = & n + \ M *" V + v ^ 2r) - M if r ' 2 - (42) 

where p — (pi — P2)/V^ and r — {f\ — r 2 )/\/2 (we have also assumed a local two-body interaction). 
By solving the two-body Schrodinger equation using the previous Hamiltonian with high accuracy in a 
large model space, one can construct by means of Eq. (36) the approximate operator Cj^ and explicitly 
an effective interaction in a chosen model space. We note that the use of the HO piece is essential. 
Thus, it not only provides a term that depends on the particle number, but also ensures that all the 
eigenstates of Hamiltonian (42) are bound, and hence the HO basis is appropriate for the eigenvector 
expansion. Finally, it is important to note that, at the two-body cluster level, the unitary transformation 
accommodates mostly the short-range correlations, while the long-range any many-body correlations 
are included by increasing the model space or the cluster level. 

By construction, there are two limits in which the exact solutions are recovered: 

(i) working in the same cluster approximation, but increasing the size of the model space approaching 
the full Hilbert space; 
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Figure 1: The ground-state energy of the triton as a function of the shell model truncation parameter 
N max using an effective interaction derived in the two-body cluster approximation. We present several 
HO frequencies: 19 MeV (open diamonds), 22 MeV (full up triangles), 24 MeV (circles), 26 MeV (full 
down triangles), 28 MeV (open squares), 30 MeV (full diamonds), and 32 MeV (open triangles), respec- 
tively. The dashed line is the exact ground state for the CD Bonn interaction, calculated independently 
in a Fadeev approach [3]. Figure taken from Ref. [39], courtesy P. Navratil. 

(ii) keeping the size of the model space constant, but increasing a up to A. 

In practical applications, the cluster approximation was never applied for a > 3, and one must rely on 
property (i) in order to obtain solutions that do not depend upon parameters of the calculation. Thus, 
the dependence on the HO frequency and model space, introduced by the cluster approximation, is 
minimized by increasing the model space. Because many-body terms in the interaction are neglected in 
this approximation, one looses the variational character of the calculation and it is not unusual to have 
a convergence to the exact eigenenergy from below. However, the advantage is that one usually obtains 
a much faster convergence than if one uses the "bare" interaction, which does preserve the variational 
character. 

In Fig. 1 we illustrate the convergence properties of the effective interaction in the two-body cluster 
approximation for 3 H. In this case, the starting Hamiltonian contains, in addition to the kinetic energy, 
only the non-local NN potential CD Bonn [3]. In small model spaces, the energy can deviate quite 
significantly from the exact solution obtained in a charge-dependent Fadeev calculation in 34 channels 
[3]. However, increasing the model space to large enough values, the ground-state energy converges to 
the Fadeev approach (note that because of missing three-body interactions, the calculated value misses 
the experiment). The reason is that by increasing the model space, the neglected contributions from 
higher order clusters become smaller and smaller. In fact, it was demonstrated analytically that, if 
one starts with a "bare" NN interaction, in the large model space limit, the two-body cluster effective 
interaction is exact to second order in perturbation theory [49]. Additionally, one often can find a 
special frequency for which the rate of convergence is very large for some states, as shown in Fig. 1 for 
oo = 26 MeV. In this case, the contribution from neglected three-body clusters is small. This usually 
happens when the length associated with the HO frequency is close to the size of the respective state. 

In order to demonstrate the effect of the higher cluster approximation, in Fig. 2 we compare the 
convergence rates of the energies of the two lowest spin and isospin zero states in 4 He, obtained using 
effective interactions calculated in two- and three-body cluster approximations respectively. Like in 
the previous example, the same charge-dependent CD Bonn potential was employed. As expected, a 
higher-body cluster approximation includes more correlations in the interactions, and the convergence is 
faster. This is illustrated in the right panel of Fig. 2, where we plot the ground-state energy dependence 
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Figure 2: The ground- and first excited-state energies of 4 He as a function of the model space truncation 
parameter N max in the three-body cluster approximation (left panel), and comparison of the convergence 
rates of the ground-state energy for the two- and three-body cluster approximations (right panel) for two 
different HO energies . The dashed line is the exact ground-state energy for the CD Bonn potential used 
in this investigation, while the dotted lines represent the experimental ground- and first excited-state 
energies. Figure from Ref. [50]. 



on N max obtained by computing the effective interaction using both the two- and three-body cluster 
approximations. The rate of convergence is faster in the three-body cluster approximation for both HO 
energies chosen for this example. Like for the 3 H energy shown in Fig. 1, HO frequencies for which 
the HO length is of the size of the states considered show faster convergence to the exact solution. The 
latter differs from the experimental value mostly because we neglect three-body interactions. Unlike 
for the ground state, the first + excited-state energy has a faster convergence rate for ou — 19 MeV. 
However, this state has a much more slower convergence rate than the ground state, and even in the 
largest model spaces the results are quite sensitive to the choice of the HO energy parameter. 

The NCSM approach described here has been applied successfully to the description of a large 
number of properties of low-lying states in light nuclei [12, 51, 52, 53, 54, 55, 56, 57], nuclear reaction 
observables [58, 59] and even to predict properties relevant to the physics beyond the Standard Model 
[60]. We will review select applications below. 

The cluster approximation employed in these calculations is not justified a priori. Mintkevich and 
Barnea have demonstrated that the higher order clusters are negligible in the limit of the model space 
approaching the "full" Hilbert space [49], as expected. They have also presented arguments suggesting 
that, at least for a system of bosons, the induced many-body terms fall much faster than the effective 
two-body interactions as the model space is increased to the full space. Nevertheless, there is no clear 
limit when the model space is large enough to ensure that the neglected many-body terms are small. 
All the results, however, suggest that, as long as the procedure is used to converge with increasing the 
size of the model space, one can obtain the exact result for a given Hamiltonian, as demonstrated in a 
benchmark calculation of the 4 He ground-state properties by several few-body methods [14]. Sometimes 
one can rely on extrapolation methods to inferrer the converged value [61], although it might not be 
obvious how large the model space has to be in order to predict the correct value. 

Several other methods based on unitary transformations have been proposed in the last decade to 
derive interactions better suitable for many-body calculations [62, 63, 64], either by means of exact 
[65, 66, 67] or approximate [68, 69] methods. In this case, the goal is to "soften" the interactions, so 
that even initially the model and excluded spaces are not as strongly coupled as when using "bare" 
interactions. The similarity renormalization group (SRG) in particular provides an elegant way to 
accomplish the decoupling of the low- and high-momentum spaces, and still preserve the accuracy of the 
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initial interactions. But such approaches suffer from the same shortcoming as the unitary transformation 
in NCSM: the procedure induces up to A-body interactions and out of practical considerations one needs 
to truncate the interactions to a level where calculations can be performed, hitherto by including up to 
three-body induced (and genuine) interactions. However, as in the case of the cluster approximation for 
the unitary transformation in NCSM, there is no proof that the neglected terms are small necessarily. 
For light nuclei has been found that the contribution of four- and higher-body induced interactions 
is negligible [65, 66, 67], but calculations of medium-mass nuclei have suggested that, although the 
hierarchy of the many-body forces seems to be preserved, the neglected terms become increasingly 
important [66, 67, 70]. Finally, it is important to note that in the calculations involving SRG and 
similar approaches, one usually uses inter- nucleon interactions that are derived from EFT [71, 72]. 
Such interactions have intrinsic systematic errors that come from the neglected terms at the order 
of the expansion. And while one can argue that the similarity/unitary transformations preserve the 
accuracy of the initial EFT Hamiltonian, they also induce terms not included in the EFT expansion, 
thus potentially enhancing those contributions. Hence, the net result can be an enhancement of the 
terms that are small (to some order) and, possibly, a loss of predictability. 



3.3 Effective operators in NCSM 

In a consistent approach, the wave functions obtained using a transformed Hamiltonian could be used 
to compute other observables only if the operators associated with those observables are transformed 
similarly to the Hamiltonian. In Sec. 3.1 we derived a unitary transformation that depends only on 
the Hamiltonian, ensuring that the effective operators are energy independent. The advantage of the 
unitary transformation approach is that any general operator preserves its symmetry properties. 

Of extreme importance are the electromagnetic properties, which are well known experimentally. 
In the phenomenological shell model, a long-standing problem is the relatively large effective charges 
that were found essential in the overall description of the transition strengths. They arise from the 
truncation of the model space and previous attempts at describing them in perturbation theory have 
proved unsuccessful [73]. In NCSM, Navratil et. al have reported effective charges consistent with 
the phenomenological values [74] . A fundamental derivation within the NCSM of the phenomenological 
effective charges represents, in addition to consistency, a strong motivation to pursue the renormalization 
of general operators. 

Because the implementation of the unitary transformation for general operators is a non-trivial task 
[75, 76, 77], for non-scalar operators it was done somewhat later than the introduction of the unitary 
transformation in NCSM [47, 48]. The difficulty stems from the fact that, for a general non-scalar 
operator of rank AJ and AT, 0^ AJ ' AT \ the calculation of effective operator matrix elements require in 
general the use of transformations in different channels, i.e., Eq. (34) becomes 



^(AJAT) = PjT + &JT q(AJAT) P J>T> + &J>T' , , ,, , 



\J PjT + &\t&JT \J P'j'T' + &ji T /&J>T' 



where J, T, J', T' fulfill the usual angular momentum rules \ J — J'\ < A J and |T — T'\ < AT. We 
emphasize that the transformation operator Cj jt is determined by the decoupling conditions imposed 
on the Hamiltonian. 

As in the case of effective interactions, the effective operators will acquire A-body contributions, 
even if the original operator is one-body. However, given that the transformation can be calculated only 
approximately, higher order terms are neglected. In particular, given the complexity of the procedure, 
the effective operators have been implemented only in the two-body cluster approximation [75, 76, 77]. 

For the one-body operators special care has to be taken in order to eliminate higher order contribu- 
tions to the effective operators. Thus, in the two-body cluster, the unitary transformation writes 
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X 2 ^exp(-G 2 ), (44) 

where G 2 = (arctan(u)jj — u}\j)j with calculated at the two-body level. Using the operator 
identity 

exp(-G 2 )0 exp(G 2 ) = + [0, G 2 ] + 1 [[O, G 2 ], G 2 } + ... (45) 

one can obtain immediately for a general one-body operator = J^i^i, the expression of the effective 
operator in the two-body cluster approximation 

i i,j 

where we have retained only one- and two-body terms in the expression. 

Similarly one obtains the expression for a general two-body operator 0^> = £V . 

P 2 Of n P 2 = P 2 J2 e- Gij O ije G -P 2 . (47) 

Combining Eqs. (46) and (47) for the special case of a one- plus two-body Hamiltonian, we recover the 
expression of the effective Hamiltonian (39) in the two-body cluster approximation. 



3.4 Selected results 

NCSM has been successful in describing a large number of nuclear properties, from energy levels to 
reaction observables. While the main scope of the paper is to discuss the effective interactions used 
within the NCSM framework, we find appropriate to illustrate with a few applications to nuclear systems 
how well this approach works, especially in the case of nuclear spectra. The interested reader can find 
more examples in a recent review of the latest NCSM calculations [13]. 

We have already discussed the two- and three-body cluster approximations for the effective interac- 
tions in the case of triton and alpha particle. However, because the main goal was a demonstration of 
the efficacy of the effective interactions in numerically accessible model spaces, as well as a comparison 
of the two- and three-body approximations, the starting Hamiltonian included two-body terms only. 
As a result, we have seen that the ground-state energies do not converge to the experimental values 
in either of the cases. In addition, precise calculations have shown that two-body forces only cannot 
account for the experimental ordering of some states in a number of p-shell nuclei [78] and in general 
predict binding energies that are too small [2, 48, 79]. Today it is generally accepted that, in order 
to obtain a good agreement with the experimental data it is important to include three-body forces in 
the calculations. Experience with phenomenological forces [2], as well as rigorous analytical proof in 
the case of EFT interactions [80], has shown that four- and higher-body forces are highly suppressed, 
although it is conceivable that for some processes involving large enough momenta they could become 
significant. 

Early light nuclei calculations have been based on phenomenological two-body forces, which provide 
an excellent description of the NN scattering data up the pion threshold, augmented by three-body 
interactions. In general, the phenomenological two- and three-body interactions provide very good 
agreement for a large number of low-energy observables [2], but in this approach the relation to the 
theory of strong interactions, QCD, is completely lost. The EFT interactions, and in particular chiral 
interactions, have emerged with the potential to provide a unified approach to the derivation of two- and 
many-body interactions. The symmetries of QCD, such as the approximate chiral symmetry of QCD, 
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play a central role in the derivation. In the EFT approach, one starts with the most general Hamiltonian 
with the appropriate degrees of freedom and symmetries [4, 5]. Interactions among nucleons consist of 
pion exchanges and contact interactions, which model short-range dynamics (e.g., exchange of heavier 
mesons) and at very low energies even pion exchange can be treated as short ranged, so that the theory 
contains only contact interactions. Because the number of terms allowed by the symmetry is infinite, 
it is important to order these terms so that observables can be calculated in an expansion in powers of 
a small parameter, given by the ratio of the relevant momentum Q divided by Mqcd ~ 1 GeV. This is 
called power counting. A truncation of this expansion at any given order must respect renormalizat ion- 
group (RG) invariance except for small errors contained in higher orders. In nuclear physics, the leading 
order (LO) terms must contain non-perturbative physics to obtain nuclear bound states and resonances, 
while subleading-order corrections, if truly corrections, should be treated in perturbation theory. The 
interested reader can find more details in other publications (see, e.g., [5]). 

In most present applications one uses the so-called "Weinberg power counting," in which all the irre- 
ducible diagrams are summed to infinite order. Although this approach is disputed [81], it is practical. 
Furthermore, there is indication that for small enough cutoffs (about 500 MeV), necessary for the renor- 
malization of the Lipmann-Schwinger equation, the shortcomings of this approach can be minimized 
[82]. In Ref. [53], Nogga et. al. have used the two-body Idaho-N3LO interaction [71], derived in the 
"Weinberg" scheme, as an input to NCSM in order to compute low-energy properties of 7 Li. In addition 
to the two-body forces, with contributions up to (Q/ Mqcd) 4 , the three-body interactions appearing at 
the (Q/M Q cd) 3 order [72] have been included. The later contain, in addition to the two pion-exchange 
term, two three-body contact interactions with or without one-pion exchange involving two undeter- 
mined low-energy constants, cd and ce respectively, that need to be adjusted to experimental data. 
Two determinations have been done for these strengths: (A) a fit to the triton binding energy and and 
the N-deuteron doublet scattering length and (B) a fit to the 3 H and 4 He binding energies. Three sets 
of calculations have been performed in NCSM: one in which only two-body forces have been included 
and two in which, in addition to the two-body interactions, the two determinations for the low-energy 
constants (A) and (B) in front of the contact terms have been used. The effective interactions derived 
within the three-body cluster approximation described in Sec. 3.1 were then used to predict the low- 
lying states in 7 Li. The best agreement of the predicted binding energy with experiment has been found 
for the case (A) (about 38 MeV). As shown in Fig. 3, all three Hamiltonians used as input predict 
the correct ordering of the states. However, some splittings between low-lying states are significantly 
influenced by the choice of the starting Hamiltonian. In particular, the splitting between 7/2~ and 5/2 _ 
is best described by choice (B), in which the binding energy is only 36.7 MeV (experimentally, the 7 Li 
binding energy is 39.2 MeV). While these results might be viewed as a undesirable model dependence, 
we point out that small deviations are expected due to the exclusion of higher order terms and should 
be reduced even further when higher order are included in the initial interactions. 

The freedom of choice regarding the fitting procedure of the low-energy constants has been investi- 
gated further in Ref. [54], where cd and eg have been determined so that the binding energies of the 
A = 3 systems ( 3 H and 3 He) reproduce the experimental values. In Fig. 4 we present the prediction 
of the excitation spectra in n B for a particular choice of the combination of the two low-energy con- 
stants that also reproduces the binding energy for the alpha particle. We observe the good convergence 
properties with the model space truncation, especially for the low-lying states. This figure also illus- 
trates the overall improvement in the theoretical spectra when three-body interactions are included. 
In particular note the lowest l/2~ and 3/2" states, whose converged values are basically degenerated 
when only two-body interactions are included, while the addition of the three-body terms produces a 
splitting that is consistent with the experimental data. The same holds for the (3/2~, 5/2~) doublet. 
In addition, a large number of other energy levels and electromagnetic properties of select p-shell nuclei 
have been investigated and found in reasonable agreement with experiment [54]. 

As expected, the inclusion of three-body interactions has an important effect for the description 
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Figure 3: Excitation energy of the lowest states of 7 Li using two-body interactions only as well as two- 
plus three- nucleon interactions (3NF), with two determinations A and B of the low-energy constants 
ce and cd, compared with the experimental data. Figure taken from Ref. [53], courtesy A. Nogga. 
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Figure 4: Theoretical and experimental natural parity excitation spectra of n B. Results with NN and 
NN+NNN interactions are presented as a function of the model space truncation N max for the HO 
frequency of 15 MeV. Figure taken from Ref. [54], courtesy P. Navratil. 
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Figure 5: B(M1), in /xfy, as a function of model-space truncation N max , for the the + — > 1 + 1 transition 
in 12 C. Figure taken from Ref. [51], courtesy A.C. Hayes. 

of other observables. This was nicely demonstrated in Ref. [51], where Hayes et. al. compared the 
convergence pattern of inelastic electron scattering to the 15.11 MeV state of 12 C, muon capture to 
the ground state of 12 B, and neutrino scattering to the ground state of 12 N, explicitly showing how the 
theoretical predictions improve with the inclusion of three-body forces. Thus, as shown in Fig. 5, if 
only two-body forces are included in the Hamiltonian, the Ml strength converges to about one third the 
experimental value. However, the introduction of the phenomenological Tucson-Melbourne three-body 
interaction (TM'(99)) [83] significantly increases the strength, although the results are far from being 
convergent in model spaces accessible at the time (note the large model space dependence in this case). 
The improvement is associated with the improved strength of the spin-orbit splitting when genuine 
three-body interactions are included. 

Except for energy, most of the calculations (including those in Fig. 5) involving observables did 
not take into account the renormalization of general operators, as it is required for consistency. Before 
2004, calculations involving effective operators have been scarce and limited to scalar operators [12, 14], 
which have the same transformation properties as the Hamiltonian. The first implementation of the 
effective operator approach discussed in Sec. 3.3 to general non-scalar operators has been introduced 
in Ref. [76]. Like for the Hamiltonian, the renormalization was implemented in relative coordinates, 
as long as the operators could be written in relative coordinates [76]. As a first test, the quadrupole 
moment of the deuteron has been calculated within the two-body approximation, which is exact for 
the two-body system. Calculations in a 4ftw model space yielded 0.179 e fm 2 when the bare operator 
was used, while the value of 0.270 e fm 2 , was obtained using the corresponding effective operator in the 
same model space. The later is in excellent agreement with independent calculation of the deuteron 
quadrupole moment for wave functions obtained with the same two-body potential. Applying the same 
method to more than three particles has yielded a different result. Thus, a very small renormalization 
has been obtained even in small model spaces for B(E2). For example, in 6 Li, if a bare E2 operator is 
used B(E2; 3 + l — > 1 + 1)=2.647 e 2 fm 4 , while the same quantity in the same model space gives 2.784 
e 2 fm 4 if the renormalization procedure is performed [76]. Both calculations have used the Argonne 
V8' NN potential. Using the CD-Bonn 2000 NN interaction, the same observable calculated using the 
bare operator in 10hw yields 10.221 e 2 fm 4 , expected to be comparable with the results obtained with 
Argonne V8'. Overall, the difference between the bare operator results in the 2hu and lOftoj model 
spaces, coupled with the small renormalization at the two-body cluster level, indicate sizable effective 
many-body effects needed to correct the 2hu B(E2) value. 
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Figure 6: Left panel: renormalization of the ground state expectation value of the relative Gaussian 
operator using realistic wave functions for 4 He, as a function of the range of the operator for AMI 
(circles), Qtvjj (squares), and 8fvjj (diamonds). Right panel: expectation values in different model spaces 
for Gaussian operators of selected ranges. Figure from Ref. [76]. 



In order to better understand the renormalization behavior in the case of general operators, a test 
case has been devised in Ref. [76], where a two-body Gaussian scalar operator of variable range ao has 
been considered. Mathematically, the operator writes 



0(fi,f 2 ) 



vr 3 / 2 a3 



exp 



(4£ 



While clearly not a realistic observable, this operator can be used to observe the behavior of the 
renormalization operator with the range. The left panel of Fig. 6 illustrates this dependence, showing 
a clear contrast between the renormalization of the operators of short range, whose values change 
by up to 90% for the shortest ranges considered, and that of the long-range operators that present 
little or no renormalization. The right-hand panel demonstrates that, in the case of a short range 
operator (ao = 0.2 fm), the renormalization produces observables that are model space independent to 
a large extent, while the results obtained with the bare operator vary significantly with the cutoff. In 
contrast, for ao = 1 fm, there is little contrast between the expectation value obtained with either the 
effective or bare operator, even though, as expected, the renormalization is more seizable in the smaller 
model space. Hence, one can immediately draw the conclusion that the unitary transformation at the 
two-body cluster level strongly (and effectively) renormalizes short-range operators, while long-range 
operators are largely unaffected by the procedure. This further reinforces the assumption that the two- 
body cluster approximation renormalizes the short-range part of the interaction, while the long-range 
and the many-body effects have to be incorporated by increasing the model space. Obviously, the 
long-range operators are going to require large spaces to converge, as they are sensitive to the correct 
description of the asymptotic of the wave function. The same behavior has been confirmed within a 
SRG renormalization of general operators in Ref. [84]. 

In very few cases, it is possible to increase the model space enough so that reasonable convergence 
can be achieved, even if the operators are not consistently renormalized. This is mostly the case of the 
A = 3 system, where in relative coordinates one can increase N max to large values. As an example, we 
illustrate in Fig. 7, the dependence of the N max for the electric polarizability a# defined by 



2a 



N^0 



\{N\DMl 
En — E 



(49) 



where a is the fine-structure constant, E is the energy of the ground-state |0), E^ is the energy of the 



19 




max max 



Figure 7: The running of the 3 He electric polarizability with N max . In the left panel we show the 
results including two-body forces only, while in the right panel we include both two- and three-body 
interactions, thus obtaining an excellent description of the ground-state energy. The polarizability 
results are in good agreement with previous calculations. Figure from Ref. [55]. 

Nth excited state, \N) (all of which are in the continuum for few-nucleon systems), and D z is the com- 
ponent of the (non-relativistic, in our case) electric-dipole operator in the z direction, which generates 
the transition between those states. Two model Hamiltonians have been used in this calculation: one 
which include only EFT derived NN interactions at next-to-next-to-next-to-next order [71], and one 
where the additional three-body interactions in the local form [85] ensure a correct description of the 
binding energy. As shown, results obtained using different frequencies eventually converge to the same 
value, even if the operator has not been transformed in the same way as the Hamiltonian. In fact, as 
Fig. 6 demonstrates, for a long range operator like the electric dipole, there is little to gain by applying 
the unitary transformation. Finally, a lower bound for the error can be estimated by using the small 
variation of the observable in the largest model space for different frequencies. 

In most of the cases, however, it is impossible to follow the same kind of convergence for all the 
observables. An assessment regarding the convergence of a certain observable is often performed by 
taking into account the change with respect the size of the model spaces. Extrapolation procedures, 
based on simple assumption of an exponential dependence on N max , have been devised for both energy 
levels [61] and observables [86]. The fitting involved in the procedure allows another estimation of the 
error for each observable considered. However, it should be noted that the error strongly depends on 
the particular form chosen for the dependence of N max , and, assuming the form to be correct, whether 
Nmax is large enough so that one can reliably extract the parameters that provide the asymptotic values. 
Furthermore, like in the cases where the error is assessed from the variation in the largest model spaces, 
this represents just a lower bound. Systematic errors that can arise, for example, from the omission 
of relevant terms in the Hamiltonian are harder to assess. In principle, EFT provides the framework 
for error estimate and systematic improvement, even though it is still difficult to estimate the errors in 
systems with more than two particles. 

3.5 Single-shell effective interactions and operators in NCSM 

The successful application of NCSM to the description of properties of light nuclei motivates the effort 
to attempt an ab initio description of nuclear properties in medium mass nuclei. In the past years, 
several many-body methods, such as the importance truncated no-core shell model [87], the ab initio 
shell model with a core [36, 37, 74], group-theory approaches [88], and others [70], have been developed 
with the goal of extending the NCSM description to medium-mass nuclei. In the following we discuss 
the NCSM with a core because the goal of the current paper is a discussion of derivation of effective 
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interactions rather than many-body truncation methods on which Refs. [70, 87] are based. In addition, 
this method is naturally connected with the phenomenological shell model discussed in Sec. 2. We 
note that, for simplicity, only two-body interactions have been included in the starting Hamiltonian. 
Hence, a direct comparison with the experimental data is a bit problematic and in order to mitigate 
this shortcoming the INOY (inside nonlocal outside Yukawa) interaction [89] has been employed. This 
NN interaction minimizes the need for three-body forces in the three-body systems, reproducing the 
NN data with the same accuracy as the other NN potentials. Even for 6 Li, INOY interaction provides 
a reasonable description of the experimental data, as illustrated in Fig. 1 of Ref. [36]. However, for a 
complete ab initio description, three-body interactions should be included in future development. 

The approach discussed in Sec. 2 has the advantage that, even for nuclei as heavy as mid pf shell 
(e.g., 56 Fe), one can perform a full space calculation with modest computing power for today standards 
and still obtain an excellent agreement with the experiment. The reason is a strong restriction of the 
model space to only one HO shell as well as a limitation of the active nucleons to a small number. The 
drawback is the phenomenological fitting of the single shell interaction used in the calculation. NCSM 
provides the opportunity to derive such an interaction from an underlying NN (and NNN) interaction, 
with the rigorous inclusion of a large amount of many-body effects (e.g., core polarization). 

While the method has been introduced by Navratil et. al. in 1996 in order to the derive electro- 
magnetic effective charges in one shell [74], a comprehensive study of the NCSM-derived single-shell 
effective interactions and operators has been performed only more recently by Lisetskiy et. al. [36, 37]. 
The method is based on a unitary transformation approach performed in two steps. The first step is 
identical to the usual NCSM approach (illustrated in the previous section), allowing the calculation 
of converged A-body states in large model spaces. In the second step, an A-bodj transformation to 
a Ohoj model space is performed, and in this model space all A-bodj correlations between nucleons 
are included. However, because most of the nucleons in Ofku will be confined to the closed shell, the 
many-body Hamiltonian will contain lower rank interactions. 

In order to better understand this method, let us take 6 Li and 6 He. In this case, after the first 
step, a six-body solution can be obtained in a large model space. Since the eigenfunctions have been 
calculated in the large mode space, using Eq. (36) one can calculate the six-body operator Cj and hence 
the effective interaction in Ofku. Like for the two-body cluster for the deuteron, the six-body effective 
interaction is exact in this case. However, in Ohu model space, for 6 Li and 6 He one has to confine two 
protons and two neutrons to the 0s shell (forming a "core"), allowing only two nucleons move in the 
p-shell, which in analogy with the phenomenological shell model, we label as "valence." Thus, the only 
interactions allowed in the model space are one- (between the valence nucleons and core) and two-body 
interactions (between the valence nucleons), even after the exact six-body cluster calculation. The same 
procedure can be applied to 5 He and 5 Li, which provide the equivalent of single-particle energies in the 
phenomenological shell model, as well as to 4 He, which provides the reference energy of the closed shell 
nucleus. 

Combining the information from A = 4, 5 and 6 nuclei, taking special care as to not double-count 
contributions, one can construct a Hamiltonian in the six-body approximation that can be further used 
in order to predict the 7 Li ground-state spectrum in OHuu. We show the ground-state energy (marked 
with square symbols) for a fixed HO frequency in Fig. 8. The convergence is displayed as a function of 
the model space used in the first transformation in order to obtain four-, five- and six-body solutions. 
The ground-state results in the conventional NCSM approach, in which the seven-body calculation is 
performed in larger and larger model spaces, is displayed with full circles connected by a full line. The 
agreement between the two methods is remarkable. Note that in this approach, the core contribution 
as well as one- and two-body terms become A-dependent, because the HO CM term added to the 
Hamiltonian, as discussed in Sec. 3, corresponds to the targeted nucleus (in this case 7 Li). In contrast, 
while in the Wildenthal interaction [32] the two-body matrix elements include a mass dependence, 
the core and the single particle energies are fixed in phenomenologically derived interactions. A similar 
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Figure 8: The ground-state energy of 7 Li calculated in the two-body valence cluster approximation 
discussed in the text (filled circles connected with full line), which include exactly up to six-body 
correlations. N max represents the truncation of the model space in which four-, five- and six-body 
solutions have been obtained. Figure from Ref. [36]. 



approach within the NCSM method, marked by full circles connected by dotted lines in Fig. 8, shows the 
need for stronger and stronger repulsive three-body forces with increasing the underlying model space, 
used to derive the effective interaction. Furthermore, although the effect is smaller, the A-dependent 
core and interactions, also provide better description of the excited state energies [36]. 

Despite the remarkable the agreement between the two theoretical predictions in Fig. 8, one can 
expect an increasing importance of the neglected many-body forces, caused by the drastic reduction 
of the model space to one single major shell. This is indeed observed for A > 7, where single-shell 
calculations suggest that three- and higher-body interactions play an increasingly important role, and 
even the partial inclusion of induced three-body forces considerably improve the description [36]. 

One of the puzzles of the phenomenological shell model is the good description of the observed 
electromagnetic transition strengths, using the full-space transition operators complemented by a simple 
rescaling of the electromagnetic charges. In order to investigate the renormalization properties of the 
quadrupole operator, the NCSM approach to the single-shell interaction was also extended to operators 
[37]. In this framework, higher-body correlations (up to six bodies, in this case) have been included in 
the procedure, like in the case of the Hamiltonian. Five-body renormalization of the 5 Li and 5 He systems 
yields the one-body contribution, while the two-body part is obtained from the renormalization of 6 Li. 
Calculations of different matrix elements in the six-body system have shown that the one-body term 
contributes significantly more than half to the total strength. In the case of isoscalar transitions, the 
two-body terms become important, thus suggesting a sensitivity to higher-body correlations in general. 
In contrast, for isovector transitions, the two-body contributions remain relatively small. 

The NCSM approach to single-shell interactions and operators could be further extended to different 
directions. First, for a consistent derivation of the interactions, it is important to include genuine three- 
body interactions in the initial Hamiltonian used to obtain the initial NCSM solutions in large model 
spaces. Second, the same approach can be extended to more challenging mass regions, like the sd shell, 
where the applicability of the traditional NCSM is very limited. In this case, it is conceivable that an 
implementation of the unitary transformation proposed by Johnson [38], tested for trapped cold atoms, 
and mentioned in Sec. 2, could improve the description of a large number of observables. However, 
considerable computational power is still required in order to obtain converged solutions for nuclei with 
one and two nucleons above the closed shell. 
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4 NCSM as an effective theory 



In the previous section we have discussed a derivation of the effective interactions in truncated model 
spaces based on the use of a carefully designed unitary transportation. The starting point was a given 
"bare" interaction, so that the results depend in general upon the model Hamiltonian chosen. In 
addition, because it is practically impossible to compute all the induced interactions, which contain up 
to A-body terms, one is forced to adopt the cluster approximation that neglects interactions of rank 
greater than a chosen order. Finally, the description of low-momentum observables requires large model 
spaces, which makes the application of the method to heavy nuclei extremely challenging. 

In EFTs, the restriction to a model space generates all the interactions allowed by the underlying 
symmetries [90]. In the traditional approach in a continuum basis, the particle momenta are limited 
within the restricted space, one can treat short- distance interactions in a derivative expansion, similar 
to the multipole expansion in classical electrodynamics. The coefficients of this expansion, called low- 
energy constants (LECs), carry information about the details of the short-range dynamics. LECs change 
with the model space in such a way that low-energy observables remain (approximatively) independent 
upon the size of the model space. 

Because the derivation of inter-nucleon interactions is based on the symmetries of the QCD La- 
grangean, the EFT approach provides a modern understanding of the nuclear forces at low energies 
[5], even in the absence of exact solutions from QCD. In particular, the appearance of light pions is 
the immediate consequence of the small and spontaneous breaking of chiral symmetry. Hence, it is not 
difficult to write the most general Hamiltonian with the appropriate degrees of freedom and symmetries 
[4], which is a generalization to systems with more than one nucleon of the Lagrangian used in chiral 
perturbation theory. Interactions among nucleons consist of pion exchanges and contact interactions, 
which model the short-range dynamics like, for example, the exchange of heavier mesons. For phenom- 
ena involving momenta much lower than the pion mass, even pion exchange can be treated as short 
ranged, leaving only contact interactions (and their derivatives) in the theory [91]. 

Including an infinite number of interactions is clearly not practical. However, they can be organized 
as an expansion of the relevant momentum scale over the nucleon mass, QjM^. Such an organization, 
called power counting, allows for a consistent truncation of the relevant terms in the nuclear Hamil- 
tonian. If desired, the precision can be systematically improved by adding higher order terms. This 
approach had mostly been applied in particle physics to systems where unitarity could be accounted 
for perturbatively. In nuclear physics, the leading order (LO) must contain non-perturbative physics to 
generate nuclear bound states and resonances, while the subleading-order corrections, should be treated 
in perturbation theory. 

If the power counting of the theory with pions is debated and still being explored [92], the power 
counting is well established for the theory in which the pions are integrated out (pionless theory) [91]. 
Thus, at low enough energies, P and higher partial waves can be neglected at low orders in the two- 
body sector. Given that the S-wave two-nucleon scattering length is much larger than the range of the 
nuclear force, set by the pion mass, one can formulate an interparticle potential as a series of contact 
interactions with an increasing number of derivatives [91, 93]. Applications of the pionless theory have 
demonstrated excellent results results for 2 H [94], 3 H [95], and even the ground state of 4 He [96]. With 
trivial modifications, this EFT has also proved useful for atomic and molecular systems with large 
scattering lengths [97]. While continuum momentum-space calculations are considerably simplified in 
this EFT, they are still quite involved beyond the four-body system. Attempts have been made to 
derive bulk properties of matter using a spatial lattice [98], but the limit of applicability of the pionless 
theory (which should break down if the momenta involved approach the pion mass) with increasing 
density is at present unknown. 

In this section, we formulate the NCSM as an effective theory. In such an approach, the interaction 
preserves the form in each model space, as dictated by the power counting. Truncation to a certain 
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order involves including all interactions up to a certain rank, hence justifying the cluster approximation. 
The effective operators describing interactions with external probes can be consistently treated within 
the same framework. 

Because the contact interactions are singular, an ultraviolet (UV) momentum cutoff A has to be 
introduced in order to solve the Schrodinger equation. This is natural in a continuum basis formulation 
of EFT, but in NCSM the cutoff has to be formulated in terms of the HO parameters. Thus, since 
(N max + 3/2)a; is the maximum energy allowed in the truncated model space, we define A as the 
momentum associated with this energy in the relative coordinate of two particles: 

A = ^M N {N max + 3/2) w. (50) 

Additionally, the HO frequency sets the spacing between HO levels and provides an infrared energy 
cutoff hw or, equivalently, a momentum cutoff A = \fM N hixi. Other authors have defined the infrared 
momentum which corresponds to the maximal radial extent needed to encompass the many-body system 
to be described [16, 67], with demonstrated benefits [99]. While it would be interesting to revisit the 
calculations of Ref. [24] using the more recently proposed infrared cutoff definition, in this paper we 
restrict ourselves to A = \JM N hxi, which was initially proposed. 

A model space is defined by the two cutoffs. In the following, we will consider two types of ap- 
plications. In one type, the system is placed on a trap, and therefore the infrared cutoff is physics 
and the solutions depend explicitly upon A. In the second type of applications, we consider untrapped 
few-nucleon systems. In this case, since A = A a/ Nmax + 3/2, the running of the observables with A 
cannot be obtained by increasing fvjj arbitrarily in a fixed- N max model space, as this procedure increases 
the infrared cutoff as well, introducing additional errors. Instead, we verify explicitly that cutoff de- 
pendences decrease with increasing A and decreasing A, and we remove the influence of the infrared 
cutoff by extrapolating to the continuum limit, where fru — > with N max — > oo so that A is fixed. 
Traditional shell-model calculations use larger values for the HO frequencies, of the order of 41 A4 1 / 3 
MeV, but in this approach we are interested in a small infrared cutoff limit, which removes the HO 
frequency dependence. 

The EFT method of deriving effective interactions can be applied to either pionfull or pionless EFTs. 
For simplicity, all the applications to nuclear systems and cold atoms have been limited to pionless EFT. 
In the following we discuss various applications to nuclear and atomic systems. 

4.1 Light nuclei in EFT framework 

Attempts to introduce EFT methods into the shell model date back to 2000 [100, 43]. Such an approach 
in based, like the unitary transformation method, on underlying NN (and NNN) interactions, thus 
suffering from the same model dependence like the conventional NCSM. While initially only a contact- 
gradient expansion modeled after EFT has been employed in the derivation of the effective interactions 
[100], the approach was later supplemented by the addition of an exact summation of the relative 
kinetic energy, which accounts for the long-range behavior coming from the excluded space. The result 
is an energy- dependent effective interaction, which, to our knowledge, has been applied only to the 
description of the deuteron and three-particle systems [43, 44]. In this section, however, we present 
an approach based entirely on EFT principles, in which the Hamiltonian is defined according to power 
counting directly in the model space, without knowledge of the excluded space whose contribution is 
encoded in the LECs. 

In LO pionless theory, the Hamiltonian is composed of the relative kinetic energy and two two-body 
contact interactions in the 3 S\ and 1 S'o NN channels [91], with corresponding parameters Cq and C$, and 
one three-body contact interaction that appears in the NNN S1/2 channel [95], with the corresponding 
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Figure 9: Running with the ultraviolet cutoff for two observables: the first (0 + ; 0) excited state in 4 He 
(lower panel) and the 6 Li ground-state energy (upper panel). The discrete points represent calculations 
for different frequencies, denoted in the legend in MeV, while the continuous lines represent a fit to a 
linear dependence of 1/A for fixed u). The dashed line represents the u — > limit. 



strength parameter Dq\ 

hi 

+C$(u, N max ) ^2 ${n- + Cq(u, N max ) ^2 - 

[i<j]0 [i<j]i 

+D (u,N max ) ^2 - rj)S(fj - f k ), (51) 

[i<j<k] 

where [i < j] s denotes pairs of particles in the S- wave NN channel of spin s and [i < j < k] triplets 
of particles in the spin- 1/2 S-wave 3N channel. The LECs depend upon the truncation N max as well 
as the HO frequency, and, in order to achieve RG they have to be adjusted in each model space so 
that one preserves the physical observables. The form of the interaction us fixed (i.e., matrix elements 
of the contact interactions in all model spaces), which represents a clear departure from the unitary 
transformation approach, in which the structure of the interaction changes from one model space to the 
other. In the continuum basis approach, the two-body parameters are fixed to two-body observables, like 
the deuteron binding energy in the 3 Si channel, and two-body scattering phaseshifts. However, because 
the NCSM basis states are built using bound states only, the direct connection with the continuum 
observables is difficult. Only C 1 (u, N max ) can be directly fixed in each model space to reproduces the 
deuteron binding energy, so alternate few-body observables have to be considered in order to determine 
Cq(uj, N max ) and D (u, N max ). Hence, the 3 H and 4 He binding energies have been used to simultaneously 
fix the remaining LECs in each model space. 

With the coupling constants thus determined, one can solve the many-body problem, predict other 
states, and compare against the experimental data. Thus, in Fig. 9, the energy of the first (J n ;T) = 
(0 + ; 0) state in 4 He (lower panel) and the ground-state energy of 6 Li are plotted vs. the UV cutoff, for 
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fixed IR cutoff. As discussed earlier, we are interested in the limit of large A and small A. Very large 
N max calculations are prohibitive even for few particles, so one needs to rely on the extrapolation of 
the results obtained at accessible N max . We assume that the energy of the state runs linearly with the 
inverse of the ultraviolet cutoff for a fixed infrared cutoff, E(A,u) = E (u) + A(u>)/A. This particular 
choice, marked by the continuous lines in Fig. 9, is motivated by the same type of LO running of the 
bound-state energy in the two-body sector in the continuum [91, 94]. The IR dependence is removed by 
taking numerically the limit u — > for a fixed A (the dashed line in Fig. 9). Following this procedure, 
the predicted value for the the first excited state in 4 He is within 10% of the experiment. This represents 
an excellent agreement which can be easily explained by the proximity of that state to the four-nucleon 
continuum threshold, well within the regime of applicability of the pionless EFT. The 6 Li binding energy 
is about 70% of experiment [24]. 

One-photon exchanges between nucleons are non-perturbative only for momenta below aM^ m 7 
MeV, where a is the fine structure constant. Because the typical momenta of the bound states is 
considerably larger than 7 MeV, the Coulomb interactions do not appear explicitly in LO. We used 
the observed 4 He binding energy as a fitting parameter; the difference between that and a Coulomb- 
corrected value is a higher-order effect. The Coulomb contribution, however, grows with the square of 
number of protons and so increases for heavier nuclei. One even expects an improvement of the 6 Li 
results when Coulomb is included explicitly, since the fit to 4 He assumes that the Coulomb repulsion 
is divided among all nucleon pairs. Hence, Coulomb effects are growing with the square of the mass 
number rather than the square of the number of protons. The effect, however, yields an underbinding 
of 2 MeV for 6 Li, well within the expected size for subleading corrections. 

Although the LO results are not as precise as the ones obtained from phenomenological potentials, 
they are consistent with QCD, with the information about QCD contained in the parameters. Moreover, 
the results are improvable order by order and adding next-to-leading order terms and beyond should 
improve the precision, as it happens in continuum calculations. Neglecting operators beyond the LO 
is expected to produce an error of about 30%, so that the 6 Li result was considered the first successful 
application of the pionless EFT to A > 4 nuclei. Unlike for the calculations in the continuum, that there 
is an additional source of errors: the assumption that for fixed u the energy runs like 1/A, although 
a softer running cannot be excluded. Thus adding an extra log(a A)/A term leaves the 4 He result 
virtually unchanged, but moves the result for the ground-state of 6 Li within 15% of the experimental 
result. However, more investigations are necessary in order to pin down the running with the ultraviolet 
cutoff. 

Before a definitive conclusion can be drawn with respect to the applicability of the pionless EFT to 
the description of 6 Li, one needs to demonstrate systematic improvement when operators beyond LO 
are included in the calculation. In the approach presented so far, i.e., a fit of LECs to few-body levels, 
the determination of coupling constants becomes quickly impracticable. It is, therefore, desirable, to 
identify a method to connect the two-body contact parameters to the scattering (or lattice) data. This 
is indeed possible if one places the interacting nucleons in a Harmonic trap, relating the bound-state 
energy levels to the scattering phaseshifts [101], very similar to the Liischer's formula on the lattice[102]. 
We devote the rest of this section to the description of such an approach and review its applications to 
atomic and nuclear systems. 

4.2 Trapped atomic systems 

The area of ultra cold atomic systems has been the subject of great experimental progress in the past 
decades. The use of magnetic fields to create Feshbach resonances in cold, trapped atomic systems 
has opened up the possibility of dialing two-atom scattering lengths to values much larger than the 
typical range of the van der Waals potential, thus creating systems whose properties are universal, 
i.e., do not depend upon the details of the interaction. A ground-breaking achievement [103] is the 
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ability to further confine just a few atoms in nearly isolated sites of optical lattices formed by laser 
beams. At low temperatures, the lattice sites may be considered as HO traps. The atomic systems with 
large scattering length display close similarities to nuclear physics systems at low energy. Hence, such 
atomic systems can be studied with techniques developed in nuclear physics and, conversely, provide 
an excellent testing ground the development of for few- and many-body methods that can be further 
applied, with little or no change, to the description of nuclear systems at low energies. In the following 
we show applications of EFT to two- and few-fermion system in a HO trap, and then further extend 
the approach to nuclear systems. 

4.2.1 Two-fermion in a trap in a EFT approach 

We consider a non-relativistic system of two spin 1/2 particles of reduced mass /i that interact with 
each other in the S-wave only, although the extension to other partial waves is straight forward. In free 
space, at sufficiently low-energy, i. e. for k much smaller than the inverse of the range of the interaction 
R, the phase shift 5o(k) is given by the effective range expansion (ERE) [104]: 

k cot 5 (k) = — — + \r 2 k 2 + . . ., (52) 
0,2 2 

with a 2 , r 2 , ... , respectively, the scattering length, effective range, and higher ERE parameters not 
shown explicitly. The ERE is an expansion in powers of kR, which, for k <C 1 / R, is a small parameter. 
Generically, the sizes of ERE parameters are set by R, for example |r 2 | ~ R, although the most 
interesting are the fine-tuned cases in which the scattering length is much larger than the range of the 
interaction. For such systems, a bound (virtually bound) state exists close to the threshold, largely 
independent upon the details of the potential. 

The ERE produces model-independent results. Thus, if one needs to describe the phaseshift with a 
certain precision, the ERE can be truncated and the system described by a finite number of parameters. 
The precision can be always improved by adding higher orders in the ERE. Potentials that generate the 
same values for this finite number of ERE parameters cannot be distinguished at the fixed precision 
level and all generate the same wavefunction for distances beyond the range of the force, r > R. 

Unfortunately, ERE cannot be used to characterize systems with more than two particles. For this, 
it is necessary to go back to the notion of potential. For particles interacting in free space, it has been 
demonstrated [91] that the ERE expansion (52) can be reproduced at each power of kR by a potential 
V constructed within EFT as a Taylor series in momentum space. In coordinate space, this expansion 
for the potential writes: 

V(r', r) = C 5(f')5(f) - C 2 { [V ' 2 5(f ')] 5(f) + S(f') [V 2 5(f}] } 

+C 4 { [V ' 4 S(f ')] 5(f) + 5(f') [V 4 5(f)j + 2 [V f2 5(r ')] [V 2 5(f)] } + ... (53) 

where C , C 2 , and C 4 are parameters, and ". . . " denote interactions that contribute at higher orders. 
The contact interactions (and their derivatives) in Eq. (53) are singular so that, when one solves the 
two-body problem with V, an ultraviolet (UV) cutoff A has to be introduced. The parameters Cj depend 
upon A in such a way that observables are UV cutoff independent, as required by RG invariance. In 
the generic situation i.e. when the size of the ERE parameters are set by R, one can simply treat the 
whole potential in perturbation theory. However, in the more interesting cases when |a 2 | ^> R, like in 
nuclear physics, the C term in Eq. (53) needs to be solved exactly, while the remaining terms can still 
be accounted for in perturbation theory [91]. (Note that for the description of resonances, the first two 
terms would have to be treated to all orders.) 

Let us assume that the same system of two particles is now confined by a HO potential characterized 
by the frequency u. In the zero-range interaction approximation, the energy eigenstates e (in units of uo) 
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of this system are related to the S-wave scattering phase shift in free space 5o(k) by the transcendental 
equation [26, 101, 105, 106] 



^TZJj--{ cot W, (54) 

where b = 1/ y/JHH is defined in terms of the reduced mass and k = y/2j]we is the relative momentum of 
the two particles. Given that the trapped system has only bound states, Eq. (54) provides a connection 
between observables in the trapped system and observables in the free space. This connection allows 
us to develop a EFT approach in which two-body parameters in the expansion (53) can be fixed by 
observables in the two-body system alone, as long as one traps the system. Theoretically, one can 
always trap a many-body system by adding a CM Hamiltonian, as described in Sec. 3. However, a lot 
of experimental effort is directed toward trapped atomic systems, and, inspired by these developments, 
we first consider fermionic systems in external HO traps. For atomic systems the trap is physics, but 
when we apply the same methods to the description of self-bound nuclear systems, we remove the trap 
and take the continuum limit. This approach will be discussed later. For now, we concentrate on the 
description of trapped spin 1/2 fermions interacting via S-wave interactions. We will consider a large 
range of two-body scattering lengths a 2 , as experimentally one can tune it over a large range of values, 
including around the Feshbach resonance. We also allow for other ERE parameters, but consider those 
set by the range of the interaction. 

Using Eq. (52) the spectrum of the two-fermion system (54) in the trap can then be written as a 
function of the ERE parameters as 

r(f-f) b ( a 2 r 2 k 2 \ , x 

U 27 - ' 1 +...) (55) 



rQ-i) 2a 2 V 2 
In quantum mechanics, the two-fermion system in the harmonic trap is described by the Hamiltonian 

H = -i-V 2 + \^r 2 + V, (56) 

where the two-body interaction V is constructed in such a way so that the spectrum of H is indeed 
given by Eq. (55). In the following, we use EFT principles to construct such an interaction at low 
energies. 

Since here we are interested in systems with large scattering length, corrections beyond the lowest 
order given by Co in Eq. (53) are treated as perturbations. For that we write the Hamiltonian (56) for 
the relative motion as 

H = if (°) + + + ..., (57) 
with the wavefunction and energy decomposed accordingly, 

|^) = |^°)) + |V;( 1 )) + ^ (2) ) + ..., (58) 

and 

E = E<® + + E^ + ... (59) 

This generic expansion is valid for all eigenvalues and eigenvectors, as long as the physics of the state is 
dominated by physics beyond the range of the interaction. The superscript ^ corresponds to the order 
in perturbation theory of the different terms. The natural basis to solve this two-body problem is the set 
of HO eigenfunctions (f) n i m (f) with the corresponding eigenenergies E n = (N + 3/2)u = (2n + I + 3/2)u 
where iV is the total number of quanta and n the radial quantum number. Since we consider only 
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interaction in the S-wave, the radial part of the solution of Eq. (56) can be expanded with the radial 
S-wave HO wavef unction 

n ( r ) = ^l%-^[L^\^~ 1/ \- r2 I^L^{r 2 /b 2 ), (60) 

where is the generalized Laguerre polynomial. In this basis, the singularity of the potential (53) can 
be removed by taking into account a finite number of shells below the highest energy shell characterized 
by n = n max . As discussed in Sec. 4.1, this corresponds to having a UV momentum cutoff, 

A = }-^/2N™ ax + 3. (61) 

with N™ ax = 2n max . Equation (61) is the equivalent of Eq. (50), but it is more general as it allows 
for different particle masses. For an eigenstate E u , the corresponding wavefunction ^ v \r) of the two- 
particle system in the trap is expanded in this finite HO basis, 

^\r)=Y J ^ ) <Pn{r), (62) 

n=0 

with the coefficients to be determined. 

The eigenenergies (59) will depend on both N max and u>, E — E(N max ,u). Since N max is arbitrary, 
the energies should not depend sensitively on N max . Although this cannot be achieved in general, it 
can for the shallow levels of interest, i.e., those which are dominated by the physics at distances r > R. 
As we show in the following, this is accomplished by allowing the to depend on both N max and 
oo, = C^ u \N max ,uj). Nevertheless, at any order a residual N max dependence introduces an error 
in the calculation of shallow levels, which should be proportional to powers of 1/A. At the end of the 
calculation we want to take N max sufficiently large, A ~ 1/-R, so that the truncation error is smaller 
than the error proportional to powers of R arizing from the truncation of Eq. (53). 

The leading-order component of the Hamiltonian H (56) is 

^ (0) = -^V 2 + ^V + Cr5(f), (63) 

with corresponding wavefunctions at LO solutions of the Schrodinger equation 

(#(0) _£(«)) ^(°)(f) = 0. (64) 

By inserting Eq. (62) into Eq. (64) and using Eq. (60) one can derive a relationship between the LO 
energy E^ of any state and the coefficient CQ(n max ,u), that is [25, 26] 

i fi ^ l! 1/2) (o) 



C {0) (n u) 27r3/2& ^ — - (n + 1) ' 

In a given model space, C^\n max , uS) can be adjusted to reproduce one energy of the two-body 
system. Although we can choose any energy of the trapped two-body system, since we are interested 
in low-energy properties, we can always fix the ground-state energy in the trap 

E^ 0) =Eo(u), (66) 

so that CQ°\n max , u) is determined from Eq. (65). In particular, at unitarity (b/a 2 = 0), if we choose 
to fix the ground state energy E (u) = uj/2, the summation over n in Eq. (65) can be performed 
analytically [107]: 

C„ ) = -^ r(lw + 3/2) . (67) 
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so that in the limit n max — > oo 

jj>C { \n max ,uj)A 7T 

2^ ' ~2' (68) 

just like in the continuum calculations [91]. This limit is independent upon the value of the scattering 
length, as illustrated in Fig. 1 of Ref . [25] . 

While the ground-state energy is thus fixed to the exact value in all model spaces, the remaining 
energies will depend upon the truncation parameter n max and HO frequency, E^ = ^(n^jO;) 
satisfy Eq. (65) and in general depend not only on u but also on n max . These energies should converge 
as n max — > oo to finite values E^^oOjOj). However, for A ~ 1/-R, the theoretical errors are dominated 
by the physics of the effective range r 2 , not included in the LO. 

Physics left out at LO can be accounted for by considering higher order terms, in increasing pertur- 
bation theory order. Thus, at NLO, we include corrections as first-order perturbations on top of the 
LO wavefunction ip(°\ The NLO correction to the two-body potential in Eq. (53) comes with an un- 
known parameter C2 that can be chosen so that a second energy level is fixed, let's say the first excited 
state. However, the NLO piece induces also a nonvanishing correction to the energy used to fix Co and 
hence requires a readjustment of C so that at NLO one simultaneously reproduces the two observables 
(energy levels). While this is a valid avenue, we use a simpler method, motivated by the fact that one 
expects the NLO correction to be small. Thus, it is more convenient to split Co into a LO piece Cq°\ 
which remains unchanged, and a NLO piece c£ 1} . The latter are treated in perturbation theory, like the 
"genuine" NLO term, and adjusted so that the level already fixed at LO does not change. According 
to Eq. (53), the NLO correction to the LO potential is 

V<U = C^8(f) - Ci 1] { [V 2 5(f)] + 2 [W(f)j • V + 2<5(r)V 2 } . (69) 

The first-order corrections to the energy, E^\ and to the wavefunction, ip^(f), are obtained in first- 
order perturbation theory 

(H^-E^)^(f) = (E^-V^)^°\f). (70) 

The requirement that two energy levels have the correct values fixes the amount of change from LO, thus 
providing two equations that determine the unknown coupling constants and in each model 
space. In the case when the lowest level Eq is already fixed to a given (experimental or theoretical) 
value, Eq 1 " 1 = 0. For non- negligible range, on can alternatively choose that at LO be fixed to a level 
of the two-body spectrum without a range, while in NLO that level is shifted to the correct position 
with range, thus requiring a non- vanishing E^K These two alternatives are the HO-basis equivalent to 
fixing Cq ^ in a free-particle basis to, respectively, a known binding energy (such as the deuteron binding 
energy) or the scattering length. In either case, if we take the first excited level Ei(u) to be reproduced 
at NLO in addition to the ground state, the two equations for the determination of Cq 1 ^^-,,, u) and 
C^inmaxjw) can be written as 

E^\n max ,uj) = Ei(u) - Ef\n max ,u), i = 0, 1. (71) 

The rest of the levels are prediction, and, as in the LO, converge to the exact values, with or without 
negligible range, when n max takes large enough values, with errors dominated by 1/A. We will illustrate 
explicitly later that the magnitude of these errors is smaller than at LO, because more physics has been 
included in the theory. 

The next-to-next-to-leading order (N 2 LO) correction to the potential is given by 

V® = C {2) 5{f) - C< 2) I [V 2 5(f)] + 2 [w(f)j • V + 25(f) V 2 } 
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+Cf ] I [V 4 <5(f)j + 4 [vV 2 5(f)j • V + 4 [vW(f)j • -VV 
+4 [W(f)j • VV 2 + 25(f) V 4 }, 



(72) 



which introduces a new four-derivative parameter C A . Additionally, similar to Eq. (69), we have added 
in (72) perturbative shifts and Cf\ which can be used to compensate the energy levels fixed at 
NLO. The three new parameters can be determined so that three energy levels (e.g., the ground-state 
energy and the first two excited levels) be fixed to the correct values. 

The correction to the energy is obtained using perturbation theory up to the second order. In 
addition to the second-order correction from (69), one has the first-order correction from in 
Eq. (72): 

= (^°)|\/^ (0) > + \ {(V> (0) |^ (1) |V (1) ) + (^ {1) l^ {1) l^ (0) )} • ( 73 ) 

Taking the lowest three levels Ei(u), i = 0, 1,2, to be fixed, the system of three equations for the 
determination of CQ 2 \n max , u), (n max ,u), and Cf\n max ,u) is 

Ef\n max , u) = Ei(u) - Ef\n max , u) - E^\n max , u), i = 0, 1, 2. (74) 

The remaining levels are, again, prediction of the theory, and they eventually converge to the cor- 
rect values. If further precision is desired, one can continue the procedure to higher orders, although 
corrections beyond second-order perturbation theory are considerably more involved. 

In order to illustrate this approach, we consider now the case of two trapped particles at unitarity, 
characterized by infinite scattering length (b/a 2 = 0) and zero effective range (r 2 /b = 0). In the 
untrapped case, this situation can be realized by considering an attractive potential given by, for 
instance, a square well with a fine-tuned depth. Asymptotically, the wave function at zero energy 
behaves as 1/r and the scattering length a 2 is then infinite. In the presence of the harmonic trap, the 
only scale is set by b, or, equivalently, u, so the solutions of Eq. (55) have to be constants. The energy 
spectrum is given by the poles of the Gamma functions in the denominators, i.e., 

e n (oo) = 2n+^, (75) 

with n > an integer. At each order we use a finite number of these energies to determine the interaction 
parameters in each model space: e at LO, e and e 1 at NLO, and e , e 1: and e 2 at N 2 LO. We also 
calculate the remaining energies, which depend on n max . Lowest excited states, from the second (n = 2) 
to the fourth (n = 4), are plotted in Fig. 10 for different values of the cutoff n max in the dimensionless 
combination Ab. All results change with Ab at LO and NLO. At N 2 LO, the second excited-state energy 
is used as input to fix the coupling constants. The predicted energy levels converge as n max increases 
to the values given in Eq. (75). In addition, Fig. 10 shows that the rate of convergence with respect to 
n ma x to the exact value increases when higher-order corrections are added in perturbation theory. 

No qualitative changes are expected with respect to the unitarity case for finite scattering length. 
Thus, in Fig. 11 we consider two particles with the interaction characterized by b/a 2 — 1 and r 2 = 0. 
Like in the unitary case, the LO is iterated to all orders, whereas the higher-order corrections to the 
potential are treated in perturbation theory. The parameters at each order are adjusted so that the 
lowest levels satisfy Eq. (55) exactly. The running of the energy for selected excited states is shown in 
Fig. 11. The exact values are now slightly lower than at unitarity, as the interaction is stronger and 
the states more tightly bound. The convergence of the energy levels to the exact values improves when 
higher corrections to the potential are added, and the difference between the truncated-space energy 
and the exact result is mitigated as more corrections are included. 
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Ab Ab Ab 

Figure 10: The second (£2), third (£3), and fourth (£4) excited-state energies, in units of u, as a function 
of the dimensionless cutoff Ab for two trapped particles at unitarity. Results at LO (diamonds), NLO 
(squares), and N 2 LO (circles) are compared with the exact values (solid lines) given by Eq. (55). The 
respective coupling parameters are fixed so that the ground-state, ground- and the first excited-state, 
and ground-, the first and the second excited-state energies are reproduced in LO, NLO, and N 2 LO, 
respectively. (Note that at N 2 LO the second excited state in the leftmost panel is used as input and 
thus constant for all values of Ab.) Figure reproduced from Ref. [26], with permission from Elsevier. 

Another case, of interest not only for nuclear physics but also for cold atoms, is the presence of a 
non-negligible range. Thus, in nuclear physics, the effective range plays an important role, while in 
atomic systems it should become relatively more important as one moves away from the resonance. A 
finite interaction range R usually generates higher ERE parameters of the same magnitude, \r%\ ~ R, 
\P 2 \ ~ R 3 , etc., even when there is fine-tuning that leads to |a.a| 3> R- As a last example, we account 
for a finite effective range r 2 , for definitives choosing parameters 6/02 = 1 and r2/b = 0.1. Regardless of 
the quantity chosen as input in LO, the existence of range introduces errors that are energy dependent 
and can only be accounted for in subleading orders. As discussed earlier, we use as LO input the 
ground-state energy e given by Eq. (55), with b/a 2 = 1, but with vanishing effective range r 2 . Hence, 
the running of is the same as in the previous example. However, at NLO, we obtain Cq 1 ^ and 
C 2 from the first two states of Eq. (55) with a 2 and r 2 non-vanishing. Subsequently, the NLO is 
different from the previous section: it accounts not only for errors of order 0(a 2 k 2 /A), due to the 
explicit model space truncation, but also for implicit ones, 0(a 2 Rk 2 ) } in the potential. The N 2 LO 




Figure 11: Same as in Fig. 10, but for b/a 2 = 1. Figure reproduced from Ref. [26], with permission 
from Elsevier. 
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Figure 12: Same as in Fig. 10, but for b/a 2 
permission from Elsevier. 



1 and r 2 /b = 0.1. Figure reproduced from Ref. [26], with 



corrections are slightly more subtle [91]. The introduction of range leaves an error that can be as big as 
0{fi^r\k As ) = 0(a 2 R 2 k A ). Errors from the explicit truncation of the model space are now 0(a 2 k 4 / 'A 2 ) 
or 0(a 2 r 2 k 4 /A), and, as in general, are smaller than errors from the truncation of the expansion once 
A > l/R. These types of errors are one order in kR or k/A from NLO, which requires for control. 
In contrast, errors from the shape parameter are only 0(a 2 P 3 k A ) = 0(a 2 R 3 k A ), two orders in kR down 
from NLO. Thus, at N 2 LO we determine Cq\ C 2 , and C4 from the lowest three levels of Eq. (55), 
still with non-vanishing a 2 and r 2 and neglecting all higher-order ERE parameters. As an illustration, 
we take r 2 /b = 0.1. Qualitatively, the only change in energies with respect to previous subsections is 
the finite jump from LO to NLO due to the effective range. 

The energy levels for three excited states for two trapped particles with the interaction characterized 
by b/a 2 = 1 and r 2 /b = 0.1 are shown in Fig. 12. The lack of effective range information in LO translates 
into an asymptotic behavior that misses the correct value. However, both the NLO and N 2 LO results 
converge as n max increases, to the values given by Eq. (55) with the range, as they should. Moreover, 
the inclusion of N 2 LO corrections speeds up convergence considerably: for a fixed value of the cutoff 
n ma x the results at N 2 LO are much closer to the exact value. 

Equation (54) has been hitherto used in one direction: to adjust the LECs so that one describes a 
set of observables (energy levels) in finite model spaces. The rest of the levels constitute prediction of 
the theory, and they approach the exact values, given again by Eq. (54), in the limit of large cutoffs. 
If we now invert Eq. (54), it is possible to determine the S-wave phaseshift: 



k n cot 5 (k n ) 



, r(3/4-e n /2) 

'r(i/4-£ 6 /2) : 



(76) 



where k n = y/2e n /b is the momentum associated with the state e n . This allows us to translate the errors 
introduced by the finite HO basis in more familiar ERE parameters. As an example, let us consider two 
interacting particles at unitarity in the absence of an effective range parameter. Thus, any deviation 
from 

k cot 5 (k) = 0, (77) 

is an error introduced by n max , and even if the initial interaction is be tuned to have zero effective range, 
the truncation introduces an effective range (and other higher ERE parameters), as it can be seen in 
Fig. 13. The main purpose of NLO (and beyond) is to reduce that effective range (and the remaining 
ERE parameters). In the example considered, the model space is defined by n max = 15. Except for 
the levels used in the fit, where also the phaseshift are correctly reproduced, e n (and the phaseshift) 
deviate from the exact value in this model space. In Fig. 13, we plot kbcot8 Q {k) as a function of k 2 b 2 
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Figure 13: Two particles at unitarity: prediction for S-wave scattering phase shifts kb cot 5o(k) as 
function of k 2 b 2 in the finite model space characterized by n max = 15. The points at LO (circles), 
NLO (squares), and N2LO (diamonds) are obtained from calculated energies via Eq. (76). The dashed 
line corresponds to a linear fit of the LO curve at small kb values, which allows an estimate of the 
effective range introduced by truncation. For comparison, we also show results with the NLO potential 
fully diagonalized (black stars), which clearly illustrates the effectiveness of the perturbation theory 
treatement. Figure reproduced from Ref. [26], with permission from Elsevier. 

(it is more natural to express k in units of l/a,2, but at unitarity b provides the only length unit). At 
LO, kcot5o(k) starts off linear in k 2 (for k 2 up to 20), and then shows higher powers of k 2 as one 
approaches the cutoff. Thus, the LO truncation generates an effective range of about 0.176 in this 
particular model space, and one could also determine further parameters of the ERE. NLO and N 2 LO 
corrections reduce the size of the ERE parameters, so that the results for the phaseshifts improve order 
by order, getting closer and closer to the horizontal axis. Since A 2 b 2 = 63, at k 2 b 2 ~ 60, the errors are 
dominated by the truncation errors, and therefore the lower orders do little to improve on the previous 
order (s). However, at low momentum, the results systematically improve, as expected. Finally, because 
in current approaches to the NN interactions one often iterates sub leading orders together with LO, we 
show also the results where NLO is fully diagonalized, finding that, indeed, the phaseshifts are more 
acurately described when NLO is treated in perturbation theory. 

In this section we have constructed an interaction within a EFT framework in the case of two 
particles trapped in a HO potential. We have shown a systematic improvement and acceleration of the 
convergence to the exact results with respect to the size of the model space as more corrections to the 
potential are added, with subleading orders treated in perturbation theory. This has been observed 
for finite and infinite scattering length CclSCS clS well clS CclSCS where the range was vanishing or taking 
a finite value, like in the case of interacting nucleons. The goal that motivated this approach, i.e., a 
method to fix the LECs only to two-body observables, has been achieved and we are thus ready to 
consider systems with more than two particles. Before applying the same method to nuclear systems, 
we consider trapped atomic systems, were a host of methods have been used to obtain few-particle 
spectra that can be used to test the accuracy of the current approach. 

4.2.2 Few-body systems 

We now consider a systems of A two-component fermions confined by a HO potential. The Hamiltonian 
that describes such a system writes: 




(78) 



i<j 
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where Hq is the intrinsic HO potential, the two-body interaction between particle i and j and 
". . ." denote three- and more-body interactions. The Hamiltonian H can be expressed in terms of the 
intrinsic Jacobi coordinates (11): 



where p& is the momentum canonically conjugated to £j. In finite model spaces, we use the EFT 
method discussed above in order to define the two-body interaction. Just as in the two-body case, 
the LO many-body Hamiltonian, containing the HO term and the LO contact two-body interaction, is 
iterated to all orders, whereas higher order corrections are treated as perturbations. RG analysis shows 
that three-body forces appear beyond N 2 LO, which is the highest order we consider here. 

Because we limit our investigations to systems of three and four particles, we use the Jacobi co- 
ordinates discussed in Sec. 3 in order to solve the Schrodinger equation. For example, in the case 
of a three-particle system we use the basis states given by Eq. (12) to expand the three-body wave 
function (we use a similar approach for four particles, although the construction of the basis is more 
involved [39]). Since the basis states (12) are eigenstates of the unperturbed Hamiltonian H , they are 
characterized by the quantum number N 3 given by 



The energy of each three-body state can be written as (A3 + 3)u, and, as discussed in Sec. 3, is used 
to define the truncated model space. In general, for a system of A particles, the truncation is done 
by introducing a cutoff N™ ax defined as the largest number of quanta in the eigenstates of H used to 
construct the A-body basis. Hence, the model space will include all the eigenstates of H with energy 
up to {N% ax + 3(A - l)/2)w. 

The truncation of the many-body space is set by N™ ax , but, as we have discussed in the previous 
section, the truncation that defines the two-body interaction in a HO basis is defined by N™ ax = 2n max 
and an exact match between the two is impossible. To some extent, we have already encountered 
this issue while discussing the cluster approximation in Sec. 3.2. In that case, there was a mismatch 
between the A-body model space and the a-body cluster model space in which the effective interaction 
is calculated. In practical NCSM applications, one chooses the truncation in such a way that the 
many-body space is the minimal required to include completely the a-body space. For example, if 
we consider just two-body interactions, in particular S-wave only, N™ ax = N™ ax when one describes 
positive-parity states, and N^ iax = N™ ax — 1 for negative-parity solutions. However, in our approach 
the renormalization of the interaction in the two-body system assumes the the states lying above the 
cutoff N™ ax have been "integrated out" rather than simply discarded. Hence, their effects are implicitly 
included in the effective two-body interaction. When these two interacting particles are embedded in a 
system with a larger number of particles, the spectators will carry energies associated with the HO levels 
they occupy. For example, of the (^3 + 3)0; total energy of one of the basis states (12), (2n 2 +/2 + 3/2)co> is 
carried by the relative motion of the spectator. As such, the maximum energy available to the two-body 
subsystem is smaller than that allowed by the A-body cutoff N™ ax and some of the states removed by 
the truncation will not be accounted for by the renormalization. One way to correct for this is to use 
the interactions renormalized with a state-dependent two-body cutoff N™ ax = N™ ax — (2n 2 + / 2 ), as first 
suggested within a conventional NCSM approach in Ref. [108]. However, the resulting state-dependent 
interaction is difficult to handle in Jacobi coordinates for systems with more than three particles, and 
cannot be incorporated in a Slater-determinant basis. In order to account for all the two-body physics 
beyond our cutoff without the use of such an interaction, we simply decouple the cutoff of the many-body 
problem from that of the subcluster defining any interaction. Such a prescription has some similarity 




(79) 



A3 = 2m + h + 2n 2 + l 2 - 



(80) 
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to the truncation used in Ref. [107]. For the three-body system, we have compared our approach of 
decoupling the two- and many-body spaces and the one involving the state-dependent interaction, but 
found no evidence that the latter is more suitable. 

Each calculation is then characterized by two cutoff parameters: N™ ax for the two-body subsystem, 
and N™ ax for the few-body system. To the order we work, no three-body forces appear (this changes 
for few-nucleon systems, where, in the pionless theory, a three-body force appears already in LO), 
so we do not need to consider a separate cutoff for renormalization of a three-body subsystem, when 
considering larger systems. Actually, for the systems considered here, only S-w&ve interactions have to 
be included up to N 2 LO. Our final results are obtained by performing a double-converging approach: 
we first increase N™ ax at fixed N™ ax until they converge, and then increase N™ ax further. For two-body 
states with N2 > N™ ax , we simply set the interaction matrix elements to zero. As we increase N™ ax 
from N™ ax = N™ ax (N™ ax = N™ ax + 1) at fixed N™ ax we observe a rapid dependence on N™ ax until 
it is somewhat larger than N™ ax , the difference reflecting the typical number of quanta carried by the 
spectators. For low-lying many-body states, further enlarging N™ ax makes little difference because in 
the many-body states added, the weight of two-body states below the interaction cutoff gets smaller 
and smaller as N™ ax increases. Having achieved results for any observable of interest that are stable 
with respect to N™ ax for each N™ ax , we can then take the limit of those values for large N™ ax . We 
illustrate this approach by presenting explicit results for the low-lying energies of trapped three- and 
four-fermion systems. Our goal is to show convergence as we increase the UV cutoff, N™ ax , as well as 
systematic improvement as the order in the EFT increases. 

We consider first the three-fermion system at unitarity (6/02 = and r2/b = 0) where semi-analytical 
results exist [109]. The ground state is characterized by total angular momentum L = 1 and negative 
parity. Figure 14 shows the convergence of the energy of this state with the truncation in the three-body 
model space, N™ ax , for two values of the two-body UV cutoff: (a) N™ ax = 10 and (b) N™ ax = 18. For 
fixed N™^, the Hamiltonian does not change, so that one expects a variational behavior of the ground- 
state energy with increasing the three-body model space. Indeed, as shown in Fig. 14, the energy 
decreases until convergence is reached for a large enough three-body model space. The value of N™ ax 
for which the convergence is obtained depends on the particular value of N™^ . Thus, for jV™ ax = 10 
the energy of the three-body ground state (in LO and corrections) does not change by more than 10~ 4 
once N™ ax > 19, while for N™ ax = 18 convergence at this level is achieved for N™ ax > 31. Even though 
for fixed the errors induced by the three-body cutoff are eliminated, the errors induced by the 

truncation in the two-body sector, where the interaction is defined, can be eliminated either by taking 
jymax ^ Q i ar g e va i ues or by adding corrections that take into account physics left out by the truncation to 
a certain order, or by combination of the two. Figure 15 shows the convergence with respect to N™ ax for 
the ground-state energy at unitarity. The LO calculation converges to the exact result [109]. Moreover, 
the result improves significantly when corrections are added to the potential: at NLO the agreement 
with the exact calculation is achieved faster than at LO, and improves still at N 2 LO. Subleading orders 
thus provide systematic improvement over LO results, with the minimal additional overhead of using 
many-body perturbation theory to different orders rather than direct diagonalization. 

The direct diagonalization of the LO Hamiltonian provides not only the ground state, but also 
excited states. Furthermore, corrections beyond leading order can be carried out similarly. In Fig. 
16 we show the running with the three-body cutoff of the energy of the first excited state with the 
same quantum numbers as the ground state, L n = 1~, for the same two values of N™ ax considered 
before, and in Fig. 17 the convergence with iV™ aa: . The same 10~ 4 precision for the same two-body UV 
cutoffs considered before is achieved at somewhat larger three-body cutoffs, N™ ax > 23 and N™ ax > 35 
respectively. Like for the ground state, for a fixed N™ ax the values of energies at all orders decrease 
until convergence is reached. Note the sharp decrease of the energy as N™ ax goes from N™ ax + 1 to 
pjmax _|_ followed by small change as the three-body cutoff is further increased. This suggests that a 
small number of quanta is carried out by the spectator, so that most of the two-body physics can be 
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Figure 14: Energy in units of the HO frequency, E/oj, of the ground state L 71 = 1~ of the A = 3 system 
at unitarity, as function of the three-body model-space size, N™ ax : (a) N™ ax = 10; (b) N™ ax = 18. 
Circles correspond to LO, squares to NLO, and diamonds to N 2 LO. The dashed line marks the exact 
value [109]. Figure from Ref. [27]. 
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Figure 15: Energy in units of the HO frequency, E/u, of the ground state L* = 1 of the A = 3 system 
at unitarity, as function of the two-body cutoff, N™ ax . Notation as in Fig. 14. Figure from Ref. [27]. 



accommodated by a relatively small three-body space. The importance of having two different cutoffs 
in the two- and many-body systems is evident in this case. Indeed, if N™ ax is fixed at N™ ax + 1 one can 
see from Fig. 16 that as corrections to the potential are added, results get worse: for both values of N 2 
the energy at NLO and N 2 LO is farther away from the exact value than the value obtained at LO. As 
Fig. 17 shows, once N™ ax is decoupled from N™ ax , the corrections to the potential again improve the 
energy systematically (except at very low two-body cutoff, i.e N™ ax < 10). Agreement with the exact 
value [109] E/tu = 4.7727 is very good: for N? ax = 22, we find at LO E/u = 4.7457, slightly below the 
value E/u = 4.8554 in Ref. [25]; at N 2 LO, E/u = 4.7721. 

The qualitative features of convergence with N™ ax and N™ ax are similar for states in other channels 
although details vary. In general, there is systematic improvement as N™ ax increases, improvement 
accelerated by the inclusion of higher-order interactions [27], as expected. 

The approach presented here is not limited to the unitarity regime. As we have seen in the two- 
body case, it can be extended to finite scattering length, as well as finite effective range. Because the 
same overall features are exhibited for finite effective range, in this paper we limit ourselves to the 
case when r 2 /b = 0; concrete examples of applications to non- vanishing effective range can be found in 
Ref. [27]. Results for the lowest L = and L = 1 states (calculated at N 2 LO) as a function of b/a 2 
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Figure 16: Same as Fig. 14, but for the first excited state with L n = 1 . Figure from Ref. [27]. 
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Figure 17: Same as Fig. 15, but for the first excited state with L n = 1 . Figure from Ref. [27]. 



are shown in Fig. 18. In the limit b/a2 — > — oo, the trapped system is non-interacting and the lowest 
state has angular momentum L — 1, corresponding to a configuration with two particles in the first 
S state, and the third in the first P level. As interaction becomes stronger, the splitting between the 
two three-body states decreases, and at b/a 2 ~ 1.5, the two levels become degenerate, and the L = 
becomes the ground state, as pointed out in Ref. [25] and confirmed in Ref. [110]. A more precise 
N 2 LO calculation puts the inversion point at b/a 2 — 1.34 [27]. Increasing even further the interaction 
strength, the ground-state energy converges to — 1/2/za 2 ,, suggesting a configuration with one particle 
moving in an S- wave around a bound state of the other two (dimer). Thus, allowing the dimer to 
form inside a wide-enough trap, b ^> a 2 , the low- lying three-body spectrum can be associated with the 
spectrum of two particles (one composite) in a trap. If the atom-dimer momenta is smaller than l/a 2 , 
Eq. (55) can be used to determine the ERE parameters for atom-dimer scattering. This method has 
been applied to the computation of the atom-dimer scattering in Refs. [27, 111], and will be used below 
to calculate the quartet neutron-deuteron scattering. 

The good agreement with semi-analytical solutions is not limited to the trapped three-body system. 
Comparable precision can be achieved for the A = 4 system. As before, we fix the value of the two-body 
cutoff N™ ax and increase the size of the many-body model space, defined here as N™ ax . We show only 
results at unitarity, where other numerical solutions exist [107, 112, 113], but finite scattering length 
and effective range can be entertained in the same framework. Figure 19 shows the convergence of 
the ground and first-excited states as a function of the two-body ultraviolet cutoff N^ 10 ^. In LO, the 
ground-state energy for _/V™ ai ' = 10 is E/u = 3.64, to be compared with E/u = 4.01 obtained in Ref. 
[25]. Corrections up to N 2 LO change the prediction to E/u = 3.52, which is in good agreement with 
previous calculations where the ground-state energy was found to be 3.6 ±0.1 [112], 3.551 ±0.009 [113], 
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G-e V = 1 




Figure 18: Energy in units of the HO frequency, E/u, of the lowest states with L* = 1 (circles) and 
L w = 0+ (squares) at N 2 LO, function of the ratio 6/0,2 between the HO length and the two-body 
scattering length, when r 2 /6 = 0. At b/a 2 — 1.34 there is an inversion of ground state. Figure from 
Ref. [27]. 
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Figure 19: Same as Fig. 15, but for the ground state (left panel) and first excited state (right panel) 
with L n = + for the A = 4 system at unitarity. Notation as in Fig. 14. Figure from Ref. [27]. 



and 3.545 ± 0.003 [107]. Like in the case of two- and three-particle systems, improvement is significant 
and systematic. 

In this section, we have considered physical systems composed of spin-1/2 fermions in external HO 
traps. In order to describe such systems, we have derived effective interactions in finite model spaces in 
a EFT framework, with corrections up to N 2 LO. To this order, the interactions are purely two-body and 
determined by the two-body scattering length a 2 and effective range r 2 . We demonstrated both good 
agreement with existing solutions at unitarity [107, 109, 112, 113], as well as systematic improvement in 
the EFT expansion. In the last part of the paper, we illustrate how the same methods can be extended 
to solving the nuclear few-body problem. 



4.3 The trapped nuclear systems and the continuum limit 

The successful application of EFT principles to the derivation of effective interactions in finite model 
spaces in the case of trapped atomic systems motivates us to pursue a similar approach for nuclear 
systems. However, there is a major difference between the two systems. In the case of atomic systems, 
the trap is physical and influences the long-range behavior of the wave function. On the other hand, 
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the nuclei are selfbound in the absence of the trap. The systems become very similar at the centre of 
the trap, where the wave function has no knowledge of the traps existence. 

The harmonic trap constitutes an essential element for the renormalization of the interaction. There- 
fore, in order to proceed with the same approach to the description of the nuclear system, one needs to 
place it in an external harmonic trap. In such a case, the spectrum of the trapped two-body nuclear 
system can be related to the nuclear phaseshifts through Eq. (54), like for the trapped atomic systems 
presented before. But because in this case the HO trap represents an auxiliary tool, one needs to take 
the continuum limit when the trap vanishes, similar to the procedure discussed in Sec. 4.1. 

Another significant difference for the nuclear case as compared to the spin- 1/2 fermion case is the role 
of three-body forces. In the absence of the trap, in order to achieve RG invariance one needs to introduce 
a three-nucleon force in the pionless EFT at LO, as discussed in Sec. 4.1. This feature is not affected 
by the presence of the trap [114], as expected from the short-distance character of renormalization. 

In this section, we present application to two- and three-nucleon systems. In order to renormalize 
the interaction, we trap the system. All the calculations are performed at finite HO frequency u, but 
the final results are obtained in the limit u — > (the only exception is the calculation of two-body 
phaseshifts, which can be obtained in the presence of the trap, assuming zero range interactions). In 
the case of many-nucleon systems, we also introduce the three-body force in the appropriate three-body 
channels. 

4.3.1 Trapped two-nucleon system 

Let us first consider two nucleons in a harmonic trap. In this case, the two-body interaction V is 
constructed as in Sec. 4.2.1 from the exact spectrum in the trap (55). For exemplification, we take the 
3 Si channel, although the treatment of the 1 So channel is similar. At LO the Schrodinger equation is 
solved exactly in a finite model space characterized by the cutoff n max and the coupling constant Cq is 
adjusted such that the ground state corresponds to the lowest energy given by Eq. (55) with only the 
scattering length in the triplet channel a-n = 5.425 fm on the right hand side. Beyond LO, corrections 
are treated as pertubations. At NLO, the coupling constants appearing at this order, C$ and C\ are 
fixed such that the two lowest energy correspond to the lowest energies given by Eq. (55) with a 2t and 
the effective r 2 t = 1.749 fm. 

As discussed in Sec. 4.2.1 once the coupling constants are fixed in each order, one can calculate 
the scattering phaseshift for a certain cutoff [26] (see Fig. 13). In Fig. 20, we plot the phaseshift for 
both the singlet and triplet configurations. At low energies, as expected, one obtains good agreement 
with the experimental data, which worsens as the energy increases, but improves systematically order 
by order, as long as the momentum of the state is well below the cutoff imposed by the model space. 

In nuclear physics, in order to achieve good accuracy it is important to include more than S waves. 
Actually, P and higher waves do come up in higher orders in pionless EFT (even though we do not 
consider those orders in this paper). But a generalization of Eq. (54) to partial waves coupled to 
arbitrary angular momentum / has been derived [115]: 



where the same conventions for b, e and k as in Eq. (54) were employed. Luu et. al. [116] have used 
the generalization (81) to extract the two-nucleon phaseshifts produced by the JISP16 NN potential 
[117], constructed via an inverse scattering approach to reproduce the low-energy NN scattering data 
with high accuracy. Figure 21 illustrate the approach for the 1 Po two-nucleon channel, where a good 
agreement between the continuum and Eq. (81) results can be observed. 

If the phaseshifts can be easily calculated in the presence of the harmonic trap, its removal is essential 
if one desires the description of the two-nucleon bound state in the triplet 3 Si channel (the deuteron). 




(81) 
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Figure 20: The two-nucleon phaseshifts as a function of relative energy in the 1 So (left panel) and 3 Si 
(right panel) channels. The ERE is marked by the continuous curves, while the dashed curves mark 
the ERE up to the effective range. The EFT results, displayed by open and filled symbols at LO and 
NLO respectively, are in reasonable agreement with ERE, especially at low energies. The calculations 
are performed at u — 1 MeV and N™ ^ = 20. Figure from Ref. [28]. 
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Figure 21: The two-nucleon phaseshifts, as function of the lab energy, for the 1 Po channel (discrete 
points) computed using Eq. (81) and compared with a continuum calculation (continuous line). Figure 
taken from Ref. [116], courtesy T. Luu. 
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Figure 22: The ground-state energy of the trapped two-nucleon system in the 3 Si channel as a function 
of the HO frequency. At LO (dashed curve) the contact parameter is fixed to the triplet scattering 
length, and thus converges to l/(m N al t ), while the range corrections introduced at NLO (continuous 
curve) bring the value very close to the experimental deuteron energy, denoted by the dotted line. 
Figure from Ref. [28]. 



The trap influences the long-range behavior of the two-body wavefunction and thus the energy of the 
state. As noted before, the trap is used to determine the LECs in each order, so all the calculations are 
performed at finite u which is then taken smaller and smaller. In Fig. 22 we illustrate the dependence 
of the lowest state of the trapped two-nucleon system as a function of the HO frequency. In the limit 
co — > 0, this state converges to the deuteron energy (marked by the dotted line), as it can be easily seen 
at NLO, while the excited states will approach a continuum spectrum. Because in LO we adjust the 
coupling constant to the triplet scattering length a 2 t, the lowest energy converges to l/(m N al t ), as in 
the continuum. 

The current approach is perhaps not the best suited to solve the two-body problem. However, the 
main motivation is the further application of the same techniques to the description of few-nucleons 
systems that we approach in the next section. 

4.3.2 The three-nucleon system 

In order to describe a system of three nucleons in a trap, we turn to the Jacobi coordinate approach 
discussed in Sec. 3, and consider first the system with total spin and isospin J n = | + and T = 1/2, 
respectively. In this case, while the three-nucleon force is not forbidden, in continuum it appears at 
higher orders than what is considered here [95], and as a consequence, the properties of the three-nucleon 
system are determined by the two-nucleon input. The same is expected to hold for the trapped system. 

The two-body interactions are constructed as described in the previous section and then proceed 
with the three-body problem. Like in the case of trapped few-body atoms, one must first fix the 
truncation of the two-body space (N™ ax ) in which one performs the renormalization, while the three- 
body model-space size defined by N™ ax is increased until the ground state energy reaches a converged 
value. Representative results for the ground-state energy as a function of the two-body model space size 
N™ ax at uj = 3 MeV are summarized in Fig. 23, while qualitative features are the same for other states 
and frequencies. Clearly the energy converges to a finite value as the two-body cutoff increases. This 
convergence is displayed in Fig. 23 and constitutes a confirmation that, as in the case of a free system, 
no three-nucleon force is needed at these orders, in this channel, to renormalize the three-nucleon system 
in the trap. 

The situation changes for three nucleons in the triton channel, characterized by total isospin T = 1/2 
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Figure 23: Ground-state energy of the trapped three-nucleon system with quantum numbers T = 1/2 , 
J 71 = 3/2 + as a function of N™ ax , for u = 3 MeV at LO (circles) and NLO (squares). Figure from Ref. 
[28]. 
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Figure 24: Ground-state energy of the trapped three-nucleon system with T = 1/2 and J n = l/2 + 
as a function of A|/mjv, for different frequencies u. Calculations are performed at LO but without a 
three-nucleon force. Figure from Ref. [28]. 

and total angular momentum J n = l/2 + . In this case, the role of the three-body forces is similar to 
that of three trapped bosons [114]. Since the renormalization concerns UV momenta, it is not expected 
to be affected by the trap. In order to check this assumption, one can perform a LO calculation in the 
trap, without the inclusion of the three-body term. As before, for a fixed two-body cutoff A2 the three- 
body model space cutoff is increased until convergence is reached. Figure 24 shows the ground-state 
energy of the trapped three-nucleon system as a function of A^/m^r, which confirms its collapse as A 2 
increases. Results for the same two-body cutoff and different values of uj are close to each other, which 
is a sign of the fact that the short-range two-nucleon interaction is much stronger than the long-range 
HO potential. Indeed, this illustrates the collapse of the three-nucleon system in this channel when 
only a two-nucleon force is included in the pionless EFT [95]. 

The strong cutoff dependence demonstrates that short- distance physics has not been accounted 
for properly, independent of the existence of the trap. The proper renormalization of the system 
can be achieved by including the non-derivative three-nucleon potential already at LO. Formally, the 
Hamiltonian is composed from the free space Hamiltonian, given by Eq. (51), and the relative harmonic 
potential term. Like in Sec. 4.1, one can choose to determine the three-body contact parameter Dq 
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so that the lowest energy of the three nucleons in the trap be fixed at the experimental value of the 
triton binding energy, E t = —8.482 MeV. Obviously, by doing so one makes an error because of the 
existence of the trap, but this error decreases when the HO frequency is taken smaller and smaller. In 
Ref. [28], the behavior of the D parameter was investigated and its behavior points to what looks like 
the beginning of the limit cycle, expected in the continuum [95]. However, the maximum cutoff for 
which the calculations were possible (A^/m^v ~ 230 MeV) was only approximately half the value where 
the second branch of the limit cycle appears in the continuum calculations [95]. 

An alternative to fixing the triton energy to the experimental value in the presence of the trap is to 
consider three-body scattering observables. We have demonstrated that the presence of the trap allows 
us to determine two-body scattering phaseshifts and hence we can use the same methods in order to 
determine the three-body parameter. In the u — > limit, the lowest state does converge to the triton. 
The excited states in the trap at small u correspond to three-body "discretized" continuum states. 
Because in free space there is no bound excited state, the lowest energy states in the trap correspond 
to the scattering of a neutron on the deuteron, which can form inside the trap. Some states correspond 
to S-w&ve scattering. This is similar to the formation of the atom-dimer system in the limit bja 2 — > 0, 
discussed in Sec. 4.2.2. Thus, the three-body parameter can be adjusted so that the nucleon-deuteron 
(nd) phaseshift extracted from the lowest excited state of the trapped system [see Eq. (76)] reproduces 
exactly the experimental data. 

In order to avoid complications with the three-body forces and the bound-state, we illustrate how 
the calculation of neutron-deuteron scattering properties proceeds in the T = 1/2, J = 3/2 + NNN 
channel. The procedure is similar to the extraction of the atom-dimer scattering parameters from 
the trapped three-body spectrum [27, 111]. Thus, the requirement that the deuteron forms inside the 
harmonic trap is that its size be much smaller than the HO length parameter, a 2 t "C b, or, equivalently, 
u <C l/(m N a2 t ) — 1.8 MeV. Because the deuteron is placed already in the trap, and thus its long 
range behavior is modified, we also consider slightly larger frequencies, up to 3 MeV. The smaller the 
frequency, the larger the errors associated with the fact that smaller UV cutoff are accessible numerically. 
At fixed u, we calculate the energy spectra E 3 . n of the three-body system, as discussed at the beginning 
of this section. Assuming the separation between the NN bound state and a third nucleon, in the limit 
of zero range interaction, one can compute the nd phaseshift proceeding in a similar manner with Eq. 
(76): 

T(3/4 - {E 3 , n - E d )/2u) 



(E 3 , n - E d ) /2u cot 5 3 ^2fi nd (E 3;n -E d )^j , (82) 



r(l/4 - (E 3 . n - E d )/2u) 

where E d is the energy of the trapped deuteron, ^, nd is the neutron-deuteron reduced mass. If the 
phaseshift 5 3 {k) is given by an ERE expansion 

k cot 5 3 (k) = -— + —k 2 + . . . (83) 
a 3 2 

in terms of nd ERE parameters 03, r 3 , etc., it is a matter of simple algebra to obtain these directly 
from the spectrum of the trapped three-body system. The two-body energy E d is obtained at LO as 
the ground state solution in the 3 S 1 channel of Eq. (55) with the ERE truncated to the first term only, 
whereas at NLO, E d is the solution of the same equation, but includes range corrections. However, we 
have to impose limitations on which states from the spectrum we can reliably use in the calculations, 
as Eq. (82) holds only in the limit of zero range interactions. The range of the nd interaction is about 
the size of the deuteron itself, lZ nd ~ a 2 t, while the relative neutron-deuteron momenta associated with 
the three-nucleon state E 3 . n is p nd = ^J2^i nd (E 3 - n — E d ). From the zero-range interaction condition, 
TZ-ndVnd <C 1, one can estimate how many states one can include in the calculation. In practice, we select 
the two lowest eigenstates and extract the value of the quartet scattering length 03^ using Eq. (83). 

The values of a 3q obtained for u = 1 MeV at LO and NLO as a function of the two-body cutoff A2 
are shown in the left panel of Fig. 25. Results for other values of the HO frequency are similar. At 
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Figure 25: Left panel: Scattering length a 3q extracted from the spectrum of the trapped three- nucleon 
system in the channel T = 1/2, J = 3/2 as function of the UV cutoff A 2 , for ui — 1 MeV: LO (circles) and 
NLO (squares). Right panel: Extrapolated values a 3(J (oo) of the quartet scattering length for different 
values of u: LO (circles) and NLO (squares). The error bars correspond to the standard fitting error. 
The horizontal dotted line marks the experimental value [118]. Figure from Ref. [28]. 



a fixed oj, the scattering length a 3q should converge as A2 is increased since the energies of the three- 
nucleon system converge. As 00 gets smaller, large values for N™ ax should be considered in order to 
reduce the errors associated with the cutoff A 2 . Since in practice there are stringent limitations in the 
values of N™ ax one can handle, an extrapolation for each value of u is needed to extract the scattering 
length a 3q (oo) which would correspond to A 2 —> 00. For this purpose the trial extrapolation function 
was proposed [28] 



1 
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a3q(oo] 



+ Af + Af 
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where pi )2 and «x,2 are parameters, which are fit to the six values of the scattering length obtained at 
the largest cutoffs. 

Results of the extrapolation are displayed in the right hand panel of Fig. (25), where a 3g (oo) is 
plotted as a function of u. For HO frequencies from about 0.4 MeV to about 2 MeV the scattering 
length is such that 7.30 fm < 03^(00) < 7.53 fm and 6.08 fm < a3^ i0 (oo) < 6.16 fm for the trial 
function (84). 

While larger traps are closer to the continuum limit, the error in the scattering length increases. 
First, as u gets smaller, the imprecision on the value of a 3q stemming from the imprecision of the 
energy (i.e., the difference between the values for finite N™ ax and N™ ax — > 00) are enhanced. This can 
be understood by noticing that in Eq. (82) the energy appears with u in the denominator. Second, 
numerical imprecision also arises as oj gets smaller since the extrapolation to 03^(00) is performed in 
these cases from data at lower A 2 . With all errors taken into account one can conclude that the results 
at NLO are in good agreement with the experimental value 03^ = 6.35 ±0.02 fm [118] and with previous 
EFT calculations [95]. 

The same approach can be extended to other processes, e.g. neutron-alpha scattering. Its application 
requires only the computation of the many-body spectra in a harmonic trap, without the difficulties 
associated with the treatment of the antisymmetrization between the scattering components. However, 
if for two particles the errors are rather small for a large range of HO frequencies, when composed 
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particles are involved it is important that the range of the interaction be smaller than the HO length, 
so that the cluster states can be properly identified. In the case of neutron-deuteron scattering it was 
important to allow for the deuteron to form inside the trap. Furthermore, the relative momenta between 
clusters has to be small enough so that Eq. (54) and its equivalents, derived under the assumption of 
zero-range interactions, remain valid. 

Finally, Eq. (83) can be used, as advertised, in order to fix the three-body contact parameter D . 
Thus, D can be adjusted so that for all model spaces, defined now by N™ ax and N™ ax , the first excited 
state of the three-body system reproduces correctly the phaseshift at the corresponding momentum. 
All the other restrictions, like the size of the HO length parameter, remain in place. Once D is thus 
adjusted, one can apply the same procedure to more complex systems. 

5 EFTs and the Gamow Shell Model for exotic systems 

Nuclei located far away from the valley of ^-stability display features that do not occur for well bound 
systems. The strong coupling to the continuum manifests in the existence of halo configurations, in 
which some nucleons orbit far away from a core of more tightly bound nucleons, and of Borromean 
systems, where removal of one nucleon is accompanied by at least one more nucleon. The neutron- 
rich Helium isotopes 6 He and 8 He are two examples of such nuclei: both are Borromean halos that 
have no bound excited states. They also exhibit the "binding-energy anomaly", i.e., higher one- and 
two-neutron emission thresholds in 8 He than in 6 He. 

A microscopic description of weakly bound/unbound nuclei requires taking into account the interplay 
between bound states, scattering states, and resonances. In other words, these systems have to be de- 
scribed as open quantum systems (OQSs), in contradistinction with well-bound nuclei, which are nearly 
isolated from the environment of scattering states and decay channels ("closed quantum systems"). A 
shell-model realization OQSs is the so-called Gamow Shell Model (GSM) [119, 120, 121, 122]. GSM 
is based on the Berggren basis [123], which consists of bound, resonant and scattering single-particle 
wave functions generated by a finite-depth potential, and it provides the mathematical foundation for 
unifying bound and resonant states - the poles of the T matrix - in the context of the Schrodinger 
equation. 

In Ref. [119], GSM has been used to describe exotic nuclei as systems of valence nucleons above a core 
using phenomenological potentials for the nucleon-core potential and schematic or phenomenological 
nucleon-nucleon interactions for the valence particles. However, even when low-momentum interactions 
obtained from realistic NN potentials are employed in order to describe the interactions among valence 
nucleons [121, 122], the approach lacks proper RG invariance. In order to address this issue, a EFT 
derivation of effective interactions in GSM model spaces, similar to the one discussed in the previous 
section has been recently proposed [124]. Furthermore, because halo configurations are characterized 
by large nuclear radii compared to the size of the tightly bound core or, equivalently, by small nucleon 
separation energy compared to the core binding energy, the physics of halo nuclei is a perfect arena 
for the application of EFT which exploit separation of scales in physical systems in order to perform 
systematic, model-independent calculations for a large range of observables [125]. By using potentials 
derived from EFT for the core-nucleon interaction and the nucleon-nucleon interaction for the valence 
particles, one do not need to assume a particular form for the interactions among constituents. Moreover, 
RG invariance is automatically satisfied, and the ultraviolet cutoff, originally introduced as an arbitrary 
separation between physics kept explicit in the theory and physics treated as short ranged, can be varied 
in a consistent manner. In the following, we briefly discuss application of halo EFT with the GSM to 
the J n = + ground state of 6 He, and for more details we refer the interested reader to the original 
publication [124]. 

Because the first excited state of 4 He is very high in energy, the J 71 " = + ground state in 6 He can be 
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described as a three-body system n + n + a. The neutrons in the halo interact with the alpha particle 
via a two-body interaction V na and with each other via a potential V nn . We denote the neutron (core) 
mass by m N (M c ) and the neutron-core reduced mass by /i Nc = mjs[M c / {m N + M c ). The potential 
between the alpha core and a neutron is constructed with EFT as described in Ref. [126]. The small 
relative momentum means that neutron and alpha particle see each other, in a first approximation, 
as elementary objects. At LO there is only one contribution, which is in the P3/2 channel, and the 
"dimeron" potential projected onto this channel can be written in the momentum space as 

k'k 

V na (k , k, ho) = — — ^ 2 , 

— * — * 

where k (k') is the incoming (outgoing) relative momentum and k = ^2fi Nc E na in terms of the total 
energy E na of the not subsystem. A and B are LECs to be determined from matching to physical 
observables. Since this interaction is singular, a regularization procedure is introduced in form of an 
ultraviolet cutoff A na . Somewhat different from the procedure followed in the previous section, this is 
here achieved by introducing a smooth regulator function 

F(x) = exp(-x), (86) 

whose role is to suppress the high-energy contributions of the potential. We thus replace the potential 
in Eq. (85) by 

V na (k', k, k ; A na ) = a(Awq) + g (Awtt )jfcg F (^ML) F (k 2 /A 2 na ) . (87) 

As discussed in Sec. 4, in order for observables to be RG invariant, the LECs A(A na ) and B(A na ) must 
depend on the ultraviolet cutoff A na . At LO, A(A na ) and B(A na ) can be fixed such that the scattering 
volume and effective "range" in the ERE in the P 3 / 2 channel exactly match the experimental values 
a na = -62.951 fm 3 and r na = -0.8819 fin" 1 [127], respectively. 

The two neutrons in the halo have sufficiently low relative momentum that meson exchange can 
be considered a short-range force. The neutron-neutron potential is thus taken from the pionless EFT 
[5, 91, 94]. At LO, the potential acts only in the 1 Sq channel and in momentum space it is simply a 
constant C . As before, the potential requires regularization, for which we continue to use the function 
F(x), but now in terms of the relative momentum between the two neutrons and an nn cutoff A nn : 

V nn (k', k; A nn ) = C(A nn ) F {k' 2 /A 2 m ) F (k 2 / A 2 nn ) . (88) 

Like in Sec. 4.3.1, the LEC C can be fixed to the neutron-neutron scattering length. 
The Hamiltonian of the n + n + a systems reads, 

2 

i=l 

where r*j is the position of neutron i — 1,2 relative to the a core, and Pi the corresponding momentum. 
This Hamiltonian is translationally invariant, the recoil term pi ■ p 2 /M c stemming from the choice of 
coordinates. The three-body equation is solved using a complete set of eigenstates of the LO potential 
V na (k ; Ana). In this particular case of the ground state of 6 He, it is sufficient to consider the set of 
continuum states along the real axis since 6 He is bound. We show in Fig 26 results of the calculations 
for the ground state of 6 He as a function of the cutoffs A nn = A na . The nearly linear dive of the ground 
state seen in Fig. 26 is reminiscent of the behavior observed with LO two-body forces in EFTs for 
systems of three bosons or three-or-more-component fermions [93]. There, the dive is even faster, more 
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Figure 26: The running of the ground-state energy of 6 He as a function of the ultraviolet cutoff A nn = 
A na . Only two-body forces have been included in the calculation. 

like quadratic in the cutoff, stemming from the strong s-wave interactions among the three pairs. In 
either case, what we see is a collapse of the ground state under short-range two-body interactions similar 
to the one first observed by Thomas [128]. It is an indication that the three-body problem has not been 
properly renormalized with only two-body interactions [93]. The cutoff dependence does not decrease 
as the cutoffs increase, as one would expected from residual cutoff dependence in a renormalized system 
that has been truncated correctly, but instead increases with positive powers of the cutoffs. The solution 
to this problem has to be found outside the two-body subsystems, which are perfectly well defined and 
well described by the EFT. Indeed, this problem disappears by adding at LO, a three-body force to 
renormalize the three-body problem [124]. 

This novel approach thus provides a unique RG invariant approach to 6 He and paves the way for 
more comprehensive studies of halo nuclei with EFT. More extensive investigation of 6 He, including 
higher-order corrections and calculation of other observables (such as the ground-state radius and the 
first excited-state energy) can be carried out. At the cost of more computational resources, other 
members of the He isotope family could be investigated as well, along the lines of Ref. [119]. More 
generally, EFT provides a consistent and systematic derivation of the interactions to be used in GSM 
calculations and should become a valuable tool in the study of other three-body resonant states, such 
as the Hoyle state in 12 C. 

6 Conclusions and outlook 

NCSM remains today, about two decades after its introduction, a powerful many-body method for 
solving the nuclear many-body problem, and not only. Its applications include nuclear spectra, electro- 
weak transitions, reaction observables, and even physics beyond standard model. Trapped atomic 
systems were also studied within its framework. 

The method is similar to the phenomenological shell model and is based on a numerical diagonal- 
ization in a finite many-body basis constructed with HO wave functions. The inherent truncation to 
finite model spaces requires the use of effective interactions, which take into account contributions from 
the excluded space, whether or not the the excluded physics is known. The main purpose of the paper 
is a review of the effective interactions and operators in finite model spaces, and not a comprehensive 
overview of all the current NCSM applications. Hence, although we have discussed selected applications 
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to energy spectra and electromagnetic observables, we have not presented more recent developments 
like the importance truncated NCSM or the implementation of resonating group methods within NCSM 
framework, with applications to nucleon-nucleus scattering or fusion observables. These are important 
developments that extend the applicability of the NCSM and will be discussed in detail in an upcoming 
review paper [129]. 

We have discussed in detail the derivation of effective interactions and non-scalar operators within 
the unitary transformation approach. The starting point is a given Hamiltonian, which includes, besides 
the kinetic term, two- and three-body interactions defined in a continuum basis. The A-hody Hilbert 
space is spit in two complementary subspaces: the model space, in which the many-body calculation 
is performed, and an excluded space, whose influence is taken into account by constructing a unitary 
transformation designed in such a way that the transformed Hamiltonian does not connect the model 
and excluded spaces. All the initial symmetries of the original Hamiltonian are preserved and, because 
of the unitarity condition, observables are not affected. The main shortcoming of this approach is that 
constructing an exact transformation is as difficult as solving the many-body problem. Therefore, the 
cluster approximation, in which the transformation is calculated for two or three interacting particles, 
and then the effective interaction used in the many-body calculation, has been introduced in practical 
applications. In this approximation, induced higher-body interactions are neglected, even though there 
is no guaranty that they are small, while contributions that initially are small or inexistent can be 
enhanced. Because the lowest approximation, i.e., the two-body cluster, integrates out mainly the short- 
range part of the interaction, short-range operators, such as the relative kinetic energy, are strongly 
renormalized, while long-range operators, e.g., quadrupole, are weakly renormalized. Implementation 
of the unitary approach to general tensor operators beyond the lowest cluster approximation becomes 
extremely challenging because the number of matrix elements explodes. All these shortcomings can 
be minimized by increasing the model space, as by construction the effective interaction approaches 
the starting interaction when the model space approaches the full Hilbert space, and so do the other 
operators. 

One of the major advantages of using NCSM is its potential to provide an ah initio description 
of heavier nuclei. Thus, the importance truncated NCSM, not discussed in detail here, can push the 
limits of applicability to medium mass nuclei. Another approach, described in detail in this paper, is 
more similar to the phenomenological shell model. However, instead of using phenomenologically fitted 
interactions, single-shell effective interactions can be derived from realistic NN and (possibly) NNN in- 
teractions in free space using a secondary unitary transformation. In addition, all the symmetries of the 
original Hamiltonian are preserved, even if part of the nucleons are frozen like in the phenomenological 
shell model. The potential of such an approach has been demonstrated in the case of p-shell nuclei, and 
could be improved by a more targeted design of the unitary transformation, as proposed for trapped 
atomic systems. However, its applicability beyond the p-shell is challenging because of the numerical 
power required to converge in the conventional NCSM approach closed shell nuclei, as well as those 
with one, two and possibly three nucleons outside closed shell. 

We have also discussed in great detail a more recent alternative to the traditional unitary transfor- 
mation approach, in which effective interactions in finite NCSM spaces are derived using principles of 
EFT. Problems such as neglecting higher order clusters in the effective interaction, or difficult renor- 
malization of general low-momentum operators can be mitigated if one formulates the problem as an 
EFT in a discrete basis. In such an approach, the form of effective interaction remains the same, set 
by the operators included in the EFT expansion to a certain order, while the strengths of the LECs 
are adjusted to reproduce experimental observables. Unlike in the conventional NCSM approach, in 
which the effective interaction converges to an underlying bare interaction with increasing the size of 
the model space, in the EFT framework it is important the behavior of observables with the variation 
of the UV cutoff. RG invariance is desired, and, at a certain order, the variation of observables with 
the UV cutoff puts a lower limit on the errors induced by neglected physics in that order. The errors 
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are reduced systematically by adding correction terms which include relevant physics or correct for 
truncation errors. 

We have illustrated the approach using an effective theory without pions, but the same general 
principles can be applied to other many-body techniques, as well as to the chiral EFT. Two alternative 
procedures of fixing the LECs have been illustrated. As a first application, we have considered the 
nuclear system in a pionless EFT approach. In LO, there are three LECs, two coupling constants in the 
NN sector and one three-body parameter in the triton channel. The strength of the contact interaction 
in the 3 Si channel has been fixed to reproduce the deuteron binding energy in all the model spaces. 
Because the deuteron is the only two-body bound state, the strength of the two-body contact term in 
the 1 So and the three-body parameters have been determined so that one simultaneously reproduces 
the triton and 4 He binding energies. Following an additional procedure of removing the infrared cutoff 
oj artifact, inherent to the HO basis, the energy of the first (0 + ; 0) state in the alpha particle has been 
determined in very good agreement with the experiment. The ground-state energy of 6 Li is predicted 
with less accuracy [24], but still within 30% of the experimental value. Despite the reasonable agreement 
for a LO calculation, in order to demonstrate systematic improvement one needs to include terms beyond 
LO. 

The lack of a larger number of bound states in the NN system makes the application of this procedure 
difficult for subleading orders, or for the EFT with explicit pions. This has motivated the search for a 
more flexible procedure, in which all the two-body parameters are fixed in the two-body system alone. 
The harmonically trapped systems present the advantage that all the states are bound, and the natural 
basis states are the HO wave functions used in NCSM calculations. Moreover, it has been shown that 
the spectrum of two particles trapped in a HO potential is entirely determined by the ERE parameters 
in the continuum [26, 101, 105, 106]. The energy levels in the trapped two-body system can be used 
to determine the LECs in a finite oscillator basis [25, 26]. In this approach, the same procedure as in 
continuum was applied: the physics at the LO was iterated in all orders, while the subleading orders 
were treated in perturbation theory. Furthermore, application to the three-body system has shown that 
the exact solution at unitarity [109] can be reproduced: the errors associated with the truncation can 
be eliminated either by increasing the cutoff [25], or by adding corrections beyond the LO [27]. The 
procedure was extended to arbitrary b/a 2 ratios, as well as small values of r 2 /b ratios or other ERE 
parameters [26, 28], of great interest for applications to the nuclear systems. 

Despite similarities between the atomic and nuclear systems, the procedure devised for trapped cold 
atoms cannot be applied directly to the nuclear many-body problem. In addition to subtleties related 
to the power counting, the nuclear system is self-bound in the absence of the trap, which modifies the 
long-range behavior. The trap plays an essential role in the renormalization procedure of the two-body 
interaction, as it connects scattering observables to bound-state energy spectra. Hence, in order to take 
advantage of the renormalization technique devised for cold atoms, the nuclear system is placed in a 
trap and all calculations are performed at finite HO frequency u. In addition, one has to take the limit 
of uj — > 0, which corresponds to the continuum. Description of the three-nucleon system properties 
has been limited to the nucleon-deuteron quartet scattering up to NLO, with results that are in good 
agreement with the experimental value and previous EFT calculations. Also, the collapse of the triton 
ground state in the presence of the trap has been investigated. In order to attack systems with more 
nucleons, the three-body force has to be included properly in the calculation. 

The approach can be extended in several directions. First, the calculations of Ref. [24] could be 
repeated in order to asses whether the reasonable agreement (in LO) with experimental 6 Li binding 
energy was accidental or could be improved at NLO, thus testing the limits of the pionless EFT with 
increasing the number of nucleons. Second, the application of the same techniques presented for the 
extraction of the quartet nucleon-deuteron scattering length could be extended to other systems in 
order to extract scattering information from bound-state physics. The nucleon-o; scattering process is 
particularly well suited to such an approach because of the lack of a bound-state system in the five- 
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body system, as well as the compact size of 4 He. Such an approach would provide an alternative to 
other methods under development [58, 59, 130, 131]. Third, the investigation of general operators is 
in order, as this constitutes a strong motivation for going to an EFT inspired approach. Finally, the 
same techniques could be extended to the EFT with pions or to other many-body methods such as the 
Gamow shell model [124] . 
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